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Abstract. We present a rigorous derivation of the Ericksen-Leslie equation starting from 
the Doi-Onsager equation. As in the fluid dynamic limit of the Boltzmann equation, we first 
make the Hilbert expansion for the solution of the Doi-Onsager equation. The existence of 
the Hilbert expansion is connected to an open question whether the energy of the Ericksen- 
Leslie equation is dissipated. We show that the energy is dissipated for the Ericksen-Leslie 
equation derived from the Doi-Onsager equation. The most difficult step is to prove a 
uniform bound for the remainder in the Hilbert expansion. This question is connected to 
the spectral stability of the linearized Doi-Onsager operator around a critical point. By 
introducing two important auxiliary operators, the detailed spectral information is obtained 
for the linearized operator around all critical points. However, these are not enough to justify 
the small Deborah number limit for the inhomogeneous Doi-Onsager equation, since the 
elastic stress in the velocity equation is also strongly singular. For this, we need to establish 
a precise lower bound for a bilinear form associated with the linearized operator. In the 
bilinear form, the interactions between the part inside the kernel and the part outside the 
kernel of the linearized operator are very complicated. We find a coordinate transform and 
introduce a five dimensional space called the Maier-Saupe space such that the interactions 
between two parts can been seen explicitly by a delicate argument of completing the square. 
However, the lower bound is very weak for the part inside the Maier-Saupe space. In order 
to apply them to the error estimates, we have to analyze the structure of the singular terms 
and introduce a suitable energy functional. 



1. Introduction 

1.1. The Doi-Onsager theory. Liquid crystals are a state of matter that have proper- 
ties between those of a conventional liquid and those of a solid crystal. One of the most 
common liquid crystal phases is the nematic. The nematic liquid crystals are composed of 
rod-like molecules with the long axes of neighboring molecules aligned approximately to one 
another. A classic model which predicts isotropic-nematic phase transition is the hard-rod 
model proposed by Onsager [T5]. Onsager introduced the notion of orientational distribution 
function and considered a mean-field model in which the rod-rod interaction was modeled 
by the excluded volume effect. Following Onsager, Maier and Saupe [TU] proposed a slightly 
modified interaction potential, now known as the Maier-Saupe potential. Doi and Edwards 
[3] extended the Onsager theory for describing the behavior of liquid crystal polymer flows. 

We use x £ Q C I 3 to denote the material point and /(x, m, t) to represent the number 
density for the number of molecules whose orientation is parallel to m at point x and time t. 
For the spatially homogeneous liquid crystal flow, the Doi-Onsager equation [4] takes 

df 1 

-^ = —n-(TZf + fTZU)-TZ-(mx K -mf), (1.1) 

where De is the Deborah number, 1Z is the rotational gradient operator (see Section 3), k 
is a constant velocity gradient, and U is the mean-field interaction potential. Onsager [18] 
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considered the potential 

U = Uf(m, t) = a I |m x m'|/(m', t)dm', 

Js 2 

where a is a parameter that measures the potential intensity. In this paper, we will use the 
Maier-Saupe potential [16] defined by 

U = a |m x m'| 2 /(m',t)(lm'. 

Js 2 

This model has a free energy 

A[f}= f (/(m,i)]n/(m,t) + J/(m, i)£/(m, i))dm (1.2) 

J§ 2 1 

as its Lyapunov functional. The chemical potential is given by 

S A 

li= — =\nf + U. 
oj 

The equation (|1.1|) can be written as 

df 1 

The stress tensor is given by 

a De = — D : (mmmm) f — (nun x TZfi)f, (1-3) 
where D = ^(k + k t ) is the symmetric part of k, and 

((•))/ = f / (-)/(m,i)dm. 

JS 2 

The homogeneous Doi-Onsager equation has been very successful in describing the prop- 
erties of liquid crystal polymers in a solvent. This model takes into account the effects of 
hydrodynamic flow, Brownian motion and intermolecular forces on the molecular orientation 
distribution. However, it does not include effects such as distortional elasticity. Therefore it 
is valid only in the limit of spatially homogeneous flows. 

The inhomogeneous flows were first studied by Marrucci and Greco [T7] , and subsequently 
by many people [HI 120] . Instead of using the distribution as the sole order parameter, they 
used a combination of the tensorial order parameter and the distribution function, and used 
the spatial gradients of the tensorial order parameter to describe the spatial variations. This 
is a departure from the original motivation that led us to the kinetic theory. Wang, E, Liu 
and Zhang [21] set up a formalism in which the interaction between molecules is treated 
more directly using the position-orientation distribution function via interaction potentials. 
They extend the free energy (jl.2j) to include the effects of nonlocal intermolecular interaction 
through an interaction potential as follows: 

A[f] = k B T [ [ /(x,m,i)(ln/(x,m,t)-l) + — l — f/(x, m, t)/(x, m, t)dmdx, (1.4) 
Jn J§2 2k B T 

where k B is the Boltzmann constant, T is the absolute temperature, and the mean-field 
interaction potential U is defined by 

U(x,m,t) = k B T I I S(x, x', m, m')/(x', m', fjdm'dx'. 

Jn Js 2 
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Here B(x, x'; m, m') is the interaction kernel between the two polymers in the configurations 
(x, m) and (x', m'). It should be symmetric with respect to the interchange of m and m', x 
and x'. B is often translation invariant and hence it can be written in the form 

B(x — x'; m, m'). 

In this paper, we take the following form as in [5j [25] : 

B(x,x,m,m)=a|mxm TT^i — F — )' 

Li 1j 

where L is the lenth of the rods, and g(x) is a radial Schwartz function with J" K3 g(x)dx = 1. 
This potential neglects the interaction between orientation and position. But is sufficient in 
many cases. The chemical potential is given by 

fi=^l = k B Tlnf(x,m,t) + U(x,m,t). 
of 

The inhomogeneous Doi-Onsager equation takes the form 
df 1 

-± + v • V/ =T^V • { (D||mm + D X {I - mm)) • (V/i)/} 
+ -^K ■ (fKfi) Ti, • (m x k m/), 

(9v 

— + v-Vv = - Vp + V-r + F e , V-v = 0. 

at 

Here Du and D± are respectively the translational diffusion coefficients parallel and normal 
to the orientation of the LCP molecule, D r = is the rotary diffusivity, V is the gradient 
operator with respect to the spatial variable x. The total stress r is the sum of the viscous 
stress r s and the elastic stress r e . There are two contributions to the viscous stress, one 
from the solvent and the other from the constraint force arising from the rigidity of the rod 
(derived in [1]), 

r s = 2?/ s D + — £ T .D : (mmmm) f, 

where D = ^(k + k t ), k = (Vv) T is the velocity gradient tensor, r] s is the solvent viscosity. 
The elastic stress r e and body force F e are given by 

r e = -(mm x 7^)/, F e = -(V/i)/. 

Let Lq be the typical size of the flow region, Vq be the typical velocity scale, To = ^ be 
a typical convective time scale. Another important time scale is the relaxational time scale 
due to orientation diffusion: T r = • The ratio of these two time scales is an important 
parameter called the Deborah number 

T r £ r Vo 

De 



T k B TL 

Let r/p = £ r , r] = r/ s + rj p , 7 = rj s /r], and Re = be the Reynolds number. We denote 

U £ (x, m, t) = / / B E (x, x', m, m')/(x', m', tjdm'dx', 
Jn Js 2 

T e e =-(mmxKn e )f, K = -(V/i £ >/, 
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where the small parameter y/e = represents the typical interaction distance and 

5 e (x,x / ,m,m / ) = a\m x m'\ 2 -^g(^—^-), 
fi £ = ln/(x,m,t) + U e (x,m,t). 

We set 

/'(x, m, t) = /(L x, m, T t), v'(x, t) = v(L x, T t)/V . 

Then the non-dimensional Doi-Onsager equation takes the following form (drop the prime 
for the simplicity): 

df e 

~M + V ' V/ = D~e V ' * ^ 7||mm + 7±(I " mm) ^ ' (V/ + fVUe) } 

+^K-(Kf + fKU £ )-K-(mXK-mf), (1.5) 

<9v T 1 — T 1 — T 

- + v ■ Vv = -V P + ^Av + • (D : ( mmmm)/ ) + ^(V ■ r| + Ff), 

where 

L De L £>e 

7 " = W D| " 7± = w^- 

The system (jl .51) has the following energy dissipation relation: 

-£(/>^+^[/]) 

+ De 2 Re ^^ E ' ( 7 H mm + ~ mm )) ' V /J e )dx, (1.6) 



where 



^[/]= / / /(x,m,t)(ln/(x,m,t)-l) + if/ E (x ) m,t)/(x,m,t)dmdx. 
We refer to [25j E3] f° r t ne numerical study and the well-posedness of the system (|1.5p . 

1.2. The Ericksen-Leslie theory. Ericksen-Leslie theory [61 [10] is an elastic continuum 
theory. The liquid crystal material is treated as a continuum and molecular details are 
entirely ignored, and this theory considers perturbations to a presumed oriented sample. 
Elastic continuum theory is a very powerful tool for modeling liquid crystal devices. 

The configuration of the liquid crystals is described by a director field n(x, t). The hydro- 
dynamic equation takes the form 

+ v-Vv = -Vp+^-Av+— -i-V-ff, (1.7) 
ot Re Re 

where the stress a is modeled by the phenomenological constitutive relation: 
Here a L is the viscous (Leslie) stress 



a = a L + a E . 



a L = a\ (nn • D)nn + a^nN + a^Nn + CK4D + a^nn ■ D + a§D • nn 



FROM THE DOI-ONSAGER EQUATION TO THE ERICKSEN-LESLIE EQUATION 



5 



with 



N=^ + v-Vn + fi-n, Q =-(k t - k). 



The six constants a±, ■ ■ ■ , a% are called the Leslie coefficients. Parodi's relation [19] gives a 
constraint for Leslie coefficients: 02 + «3 = ctQ — 05. While, a E is the elastic (Ericksen) stress 

^-W (Vn)T ' (L9) 

where Ep = Vn) is the Frank energy. The dynamic equation for the director field is 

given by 

n x (h- 71 N - 72D • n) = 0, (1.10) 

where 71 = «3 — «2 , 72 = «6 — ^5 , and h is the molecular field 

5E F „ dE F dE F 
h = — = V 



8n d(Vn) dn 

In this paper, we will consider Ep = § Jo |Vn(x)| 2 dx. In this case, we have 

h = feAn, a E = -kVnQVn = -k(Vin k Vjn k ) 3x3 . (1.11) 
The energy dissipation for Ericksen-Leslie equation is given by 
d / f Re .9 , „ 
dt\J a 2(1-7)' 



~ 2 

/ f— ^— |Vv| 2 + (ai + ^-)|D:nn| 2 + a 4 D:D 

7n V 1 - 7 71 



+ (a 5 + a 6 - — )|D • n| 2 + —In x h\ 2 )dx. (1.12) 
7i 7i y 

We refer to [5] for a derivation of (I1.12p . Concerning the mathematical study of the simplified 
Ericksen-Leslie equation, we refer to |11| IT2| IT3| [T3] and references therein. 

1.3. From the Doi-Onsager theory to the Ericksen-Leslie theory. Two kinds of the- 
ories were put forward to investigate the liquid crystalline polymers from the different points 
of view. The Ericksen-Leslie theory is phenomenological in nature, and will be become invalid 
near defects where the director cannot be defined. The Ericksen-Leslie equation contain six 
unknown parameters called the Leslie coefficients, which are difficult to determine by using 
experimental results. Especially, whether the energy defined in (|1.12p is dissipated remains 
unknown in Physics. Hence, it is very important to establish the relationship between two 
theories. 

Kuzuu and Doi [9] formally derive the Ericksen-Leslie equation from the Doi-Onsager 
equation (jl.ip . and determine the Leslie coefficients. However, the Ericksen stress is missed 
in the homogeneous case. E and Zhang [5] extend Kuzuu and Doi's formal derivation to the 
inhomogeneous case. To recover the Ericksen stress, they find that the Deborah number De 
and the interaction distance y/e should satisfy De ~ e. 

Roughly speaking, Kuzuu and Doi shows that when the Deborah number is small, the 
solution / of (jl.lj) has the formal expansion 

/ = /o(m • n) + e/H , 

where /o(m • n) denotes the equilibrium distribution function satisfying 

n ■ (nf + foRUf ) = 0, 
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and n is determined by (|2.3p . E and Zhang shows that the solution (/, v) of (jl.5p has the 
formal expansion 

/ = /o(m ■ n) + e/H , 

v = v + evi H , 

where (vo,n) is determined by (jl.7p and (jl.lOp . 

The main goal of this paper is to give a rigorous derivation of the Ericksen-Leslie equation 
from the Doi-Onsager equation. This is a singular small Deborah number limit problem. To 
justify this limit, we first make the Hilbert expansion for the solution of the Doi-Onsager 
equation, then show that the error term is small in a suitable Sobolev space. The existence of 
the Hilbert expansion is connected to the question whether the energy of the Ericksen-Leslie 
equation is dissipated. We will show that the energy is dissipated for the Ericksen-Leslie 
equation derived from the Doi-Onsager equation. The error estimates rely heavily on the 
spectral analysis of the linearized Doi-Onsager operator around the critical point, which 
includes 

1. Give a complete classification for all critical points h of A[f], which satisfies 

K-{Kh + hTUAh) = 0. 

2. The spectral analysis of the linearized Doi-Onsager operator Qh around a critical point 
h defined by 

Qhf = K ■ {TZf + hlZUf + fKUh) . 

3. Establish a precise lower bound for the bilinear form (Ghf^H^f) with 

g £ h f = n-{nf + hizu £ f + fizuh), uy = f -+u £ f. 

The first point has been given by the second author and coworkers [15]. The second point and 
the third point are completely new. To prove the second point, we introduce two important 
auxiliary operators Ah and Hh defined by 

A h f = -n-{hKf), u h = f -+uf. 

It is easy to see that Qh = —AhHh and Hh is self-adjoint. Then we reduce the spectral 
analysis of Qh to that of Tih- The proof of the third point is very subtle. Since the orthogonal 
structure is destroyed when e ^ 0, the interactions between the part inside the kernel of Q^ 
and the part outside the kernel become very complicated. To prove a lower bound, we find a 
coordinate transform and introduce a generalized kernel space of £?^(this is a five dimensional 
space called the Maier-Saupe space) such that the interactions between two parts can be seen 
explicitly by a delicate argument of completing the square. 

With the above preparations, it is still not enough to complete the error estimates in the 
inhomogeneous case. When e / 0, we can only get a strong lower bound of (Qhf, 1t e h f) for the 
part outside the Maier-Saupe space, and a weak lower bound for the part inside the Maier- 
Saupe space. In order to apply them to the error estimates, we have to analyze the nonlinear 
interactions between two parts for the singular term like ~(/r, dt(j^)fii) and introduce a 
suitable energy functional. Since we have no decay in e for the part of fn inside the kernel, 
the term dt(-^)/ii) seems to have an order i(Very singular). Surprisingly, we will 

show that it is bounded. 
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We believe that the spectral information of the linearized operator will be very important 
to study the other problems like the nonlinear stability and instability of the critical points. 
These will be left to the future work. 

2. PRESENTATION OF MAIN RESULTS 

2.1. The homogeneous case. We consider the homogeneous Doi-Onsager equation 
df e 1 

=-TZ ■ (Kf + fWUf) - K{m x (D - n) ■ m/ £ ) . (2.1) 

Here D = i(« + k t ), Q = ^(k t — k), and e is the Deborah number. The corresponding stress 
tensor a 6 is given by 

1 1 

a e = -D : (mmmm)ye (mm x IZfi 6 ) fe (2-2) 

2 S 

with n e = lnf e + Uf £ and Uf = a L 2 |m x m'| 2 /(m',t)dm'. 

In the homogeneous case, the Ericksen-Leslie equation is reduced to 

nx (-^ + n-n-AD-n) = 0, (2.3) 

together with the stress a L given by 

a L = ai(nn • D)nn + o^nN + a3Nn + Q4D + a^uu ■ D + a^D • nn. (2.4) 
Our main results are stated as follows. 

Theorem 2.1. Let h„ D be a stable critical point of A[f], and n(t) be a solution of \2. 3\) with 
the initial data no € S 2 and A given by 

<3(m-n) 2 -l) h 

A(a) = , duav — , u =Uh^ n , (2.5) 

\9o d e )hr,, n 

and go is a solution of h4-8\ ). Assume that the initial data /o(m) € -fT 1 (S 2 ) with J §2 /o( m )dm = 
1 takes the form 

3 

/ £ (m) = ^,n (m) + ^ e fc /fc(m,0) +e 2 /^ (m), 

k=l 

where fk(m,t)(k = 1,2,3) is determined by Proposition [b\l\ and/^ (m) satisfies ||/Jj IIh- 1 (§ 2 ) — 
C . Then for any T > 0, there exists an £q > such that for each < e < Eq, the solution 
/ e (m, t) of &2.1\) takes the form 

3 

f £ (m,t) = h vMt) (m) + J2e k f k (m,t)+e 2 f £ R (m,t), 

k=l 

where /jj(m, t) satisfies 

ll/M*)llif-HS 2 ) < C for any t € [0,T\. 
Remark 2.1. The Ericksen-Leslie equation \2. 3\) is equivalent to 



+ n • n - A(I - nn)D • n = 0. 

at 

Lt is easy to show that it has a unique global solution. 
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Let S 2 = (-fMm • Q ))h VI1 and S4 = {Pi(m ■ n)}^ , where Pk(x) is the k-th Legendre 
polynomial. We take the Leslie coefficients a\, ■ ■ ■ , otQ in the definition of a L as follows 

oi = -y, a 2 = -i(l + l)S 2) a 3 = -I(l-i)5 2 , (2.6) 
"4 = ^ - 2^5 2 - — S4, a 5 = -S4 + -S 2 , a 6 = -S4--S 2 - (2.7) 

Then we have 

Theorem 2.2. Let p(t) = — §D : nn. For any T > 0, i/iere exisi an £0 > suc/i that for 
each < e < £0, there holds 

\a £ (t) - a L (t) - p(t)I\ < Ce for t G [0, T]. 

2.2. The inhomogeneous case. In order to derive the Ericksen-Leslie equation with the 
Ericksen stress, we have to consider the system (|1 .5j) with De = e. For the simplicity of 
presentation, we will consider the case when the translational diffusion coefficients vanish. 
Then the non-dimensional Doi-Onsager equations takes 

^ + v £ • V/ £ = -K ■ (Kf £ + f e nU £ f) -K-(mx K £ - rnf ), (2.8) 
at e 

0v £ T 1 — "7 

— + v e • Vv £ = -Vp £ + ^Av £ + • (D e : (mmmm) /e ) 

1-7 ' 
eRe 

where k £ = (Vv £ ) T ,~D £ = \ {k £ + (k £ ) t ), and 

r £ e = -(mmx^u e ) /£ , F| = - (Vfi £ ) f e , fi £ = In f + W £ /, 



+— ^(V-r e e + Ff), (2.9) 



24/ = / / almxmf-^^^^x'.m'.^dm'dx'. 



We also require that the Fourier transform of g satisfies 

0<g(0<l for e^O, g"(0)<0. 
Now we can derive the full Ericksen-Leslie equation 

n x (h- 71 N - 72D • n) = 0, (2.10) 

+ v . V v = -Vp+-f Av + — lV-a, (2.11) 
ot ite Re 

where 71 = a 3 -a 2 , 72 = a 6 -a 5 , h = fcAn, N = |l+v- Vn+ft-n, D = ±(Vv+(Vv) T ), £1 = 
\ (Vv - (Vv) T ) , and a = a L + a E with 

a L = «i(nn • D)nn + a 2 nN + Q3N11 + CK4D + a^nn ■ D + et^D ■ nn, 

a E = -kVnQ Vn. 

Our main result is stated as follows. 

Theorem 2.3. Let Q, = M 3 and /i^n be a stable critical point of A[f], and let (n, vo) € 
C([0, T];H 20 (n)) be a solution of J£7ffl )- fO?l) on [0, T] for some T > with the initial data 
(no,vo,o), A given by \2. 5\) . and Leslie coefficients defined by /l2.6\) - (277\ ). Assume also that 
there exist constant vector c E S 2 and constant cq £ (0, 1) such that 

|n(x,t) x c\ > c for any (x, t) £ O x [0, T]. (2.12) 
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Assume that the initial data (/q (x, m), Vq) with J §2 /q(x, m)dm = 1 takes the form 

3 

/o(x,m) = ^ >no(x) (m) +^e fc / fc (x,m,0) + £ 3 /|> >0 (x, m), 

k=l 

2 

vo(x) = ^e fe v fci0 (x) + e 3 ^ i0 (x), 

fc=0 

where (/i, f2, fs, v 1; v 2 ) is determined by Proposition \7A\ and (/jj (x, m), vfj ( x )) satisfies 

WfR,o\\H 2 (nx§ 2 ) + \\ v Rfi\\H 2 (n) < C < II/r,oIIl 2 (^x§ 2 ) < Ce. 

TTien t/iere exisi £o > suc/i that for each < e < Eq, the system \2.8\) - 112. 9\) has a unique 
solution (/ e (x, m, t), v £ (x, m, t)) on [0,T] which takes the form 

3 

/ e (x, m, t) = h vMX)t) (m) + e fc /fc(x, m, t) + e 3 /^(x, m, i), 

k=l 

2 

v e (x,t) = ^e fc v fc (x,t)+e 3 v| ? (x,t), 

fc=0 

where (/j^v^) satisfies 
/or any t G [0, T] . 

Remark 2.2. T/ie non- degenerate assumption 12.12\) allows us to construct a global coordi- 
nate transformation, which is the key to establish a lower bound of a bilinear form associated 
with the linearized operator in Section 5. 

Remark 2.3. We will study the existence of the solution for the full Ericks en- Leslie equation 
in a separate paper. We refer to |12] for the simplified Ericksen-Leslie equation. 



3. Classification and stability of critical points of energy functional 
We consider the homogeneous energy functional A[f] defined by 

A[f] = [ (/(m)ln/(m) + J/(m)W/(m))dm 

JS, 2 Z 

for / G L 2 (§ 2 ). We define 

V (§ 2 ) = {<p£ L 2 (S 2 ) : J ^ <p(m)dm = o}. 

We are concerned with the local minimizer of A\f\. That is, we find all h S L 2 (§ 2 ) such 
that 

A[h + e<f>] > A[h] 

for all (j) S ^(S 2 ) when e is small enough. Taking a formal expansion, we find that 

A[h + ecf)] = A[h] + e( In h + Uh, (f) + e 2 (j- + Ucf>, <f) + 0(e 3 ). 
This motivates us to introduce the following definition. 
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Definition 3.1. We say that h £ L 2 (§ 2 ) is a critical point of the energy functional A[f] if 

— ^r-L_, = In h + Uh = const. 
df 

A critical point h is said to be stable if for any 4> G "Po(§ 2 )> there holds 

(^+U<f>,<f>)>0. 

It is easy to see that if h is a critical point of A[f], then h is a solution of stationary 
Doi-Onsager equation 

K - (Kh + hTUAh) = 0. (3.1) 

A complete classification for all critical points of A[f] was given by Liu, Zhang and Zhang 
[E]; see also [31 [71 Elj . 

Proposition 3.1. All the critical points of A[f] take the form 

e r?(m-n) 2 

^' n(m) = JJ evi^do-' 
where n is an arbitrary unit vector, and r\ = 77(a) is determined by the equation 

3 + 2r] + -!-. (3.2) 



/ e^dz « 
Furthermore, we have 

• For all a > 0, rj = 0(i.e. h = j-) is always a solution. 

• For a < a* ~ 6.731393, r/ = is the only solution. While for a = a* , there is another 
solution i] = rf . 

• For a > a*, besides rj = 0, there are exactly two solutions rj = 771(a), 7/2(0) satisfying 

- 7/1 (a) > 77* > ?? 2 (a), lhna^c,. 771(a) = lim Q ^ Q . 772(a) = rf 7 

- 771(a) is an increasing function of a, while 772(a) is a decreasing function; 

- 7/2(7.5) = 0. 

Except the monotonicity of 77(a), the others have been proved in [15]. The monotonicity 
will be proved in Lemma 14.21 

Concerning the stability of the critical point, Zhang and Zhang [22] showed 

Proposition 3.2. h = ^- is a stable critical point of A[f] if and only if a < 7.5; If a > a* , 
h m>n is stable, while h V2tIl is unstable. 

Let us conclude this section by collecting some properties of the rotational operator, which 
will be used throughout the paper. Let m £ § 2 and V m be the gradient operator on the unit 
sphere S 2 . The rotational gradient operator 1Z is defined by 

1Z = m x V m . 

Let (9, (p) be the sphere coordinate on S 2 . Then 1Z can be written as 

1Z =(— sin 4>i + cos <Aj)<9e — (cos 6 cos 4>i + cos 6 sin 0j — sin Oh) -ds 

sin 6 

d =ifti +j1Z 2 + klZ 3 . 
The following properties can be easily verified. 
1. TZ-TZ = A g2 ; 
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2. TZiiTij = —e^rrik, where m = (mi, wi2, m 3 ). If u is a constant vector, then 

1Z(m • u) = m x u, 1Z- (m x u) = —2m • u; 

3. [Kj, K k ] = j jk Ki\ 

4. jg 2 K/i/ 2 dm = - / sa /iK/sdm; 

5. = 0. 

Here A§2 is the Laplace-Beltrami operator on § 2 , e ij is the Levi-Civita symbol. 

4. Spectral analysis of the linearized operator 

We linearize the Doi-Onsager equation 1Z ■ (IZf + flZUf) = around a critical point h. 
The linearized Doi-Onsager operator Qh is given by 

g h f A =n-{nf + hnuf + fnuh). (4.i) 

We denote by H k (S 2 ) the Sobolev space on § 2 , and H^(S 2 ) = H 2 (S 2 )nV . Qh is a bounded 
operator from H 2 (S 2 ) to L 2 (S 2 ), and has the discrete spectra.This section is devoted to 
studying the kernel and spectra of the linearized operator Qh . These information will play 
a vital role in the study of small Deborah limit, and will be very important in the study of 
nonlinear stability and instability of the critical point. 

When h is a trivial critical point ho = the linearized operator Qh is reduced to 

g h0 f = A §2 (f + ±-uf). 



Aa 
5 ' 



Proposition 4.1. The eigenvalues ofQh are A^ = — k{k+l) (fork 7^ 2, k > 1) and —6 + 

and the corresponding eigenfunction is the spherical harmonics Y^/ of degree k. Specifically, 
Qh has a positive eigenvalue if and only if a > 7.5. 

Remark 4.1. The critical value 7.5 is consistent with that in Proposition 1 3.2\ deduced from 
the energy stability analysis. 

Proof. Let tp be an eigenfunction of Qh associated with the eigenvalue A, that is, 

G ho ijj = Xip. 

We choose the spherical harmonics {Y2/}i<^<5 of degree 2 as 

Yi = (m\ - m 2 ,), Y 2 = m\ - -, Y 3 = mim 2 , Y 4 = mim 3 , Y 5 = m 2 m 3 . 

Then we make a spherical harmonics expansion for ip: 

5 

ip = ^2^iYi+ Vk,iY k/ . (4.2) 

i=l k^2,l 



We have 



A §2 Uip = -aK-K / (m • m') V(m')dm' 

= —2alZ ■ / (m • m')(m x m')'i|i(m')dm' 
= 6a (m • m') 2 ^(m')dm'. 

J§2 
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Hence, 

\ip = A S 2-0 H / (m • m / ) 2 V'(m / )dm / 

2vr J S 2 



A§2'0 + —mirrijMij (4.3) 



3a 

with Mij = L 2 rrijmj'0(m)dm. Noting that 

A S 2Y k/ = -k(k + l)Y k/ , 
and plugging (|4.2p into (|4.3p . we find that 

A/ifc/ = -k(k + 

This implies that 

A = — k{k + 1) or /ifc^ = 0. 
Hence, if A > 0, then /i^ = and we have 

5 

^( m ) = ^2mYi- 

i=l 

Then it follows from ()4.3|) that 

5 

(A + 6)J2^Yi 

i=i 

= ^ / ( \y x Y{ + jW 2 ' + 2y 3 n' + + 2F 5 y 5 ' + J) V /^dm 

= p f {l^Yf + jW 2 Y 2 ' 2 + 2fi 3 Y 3 Yf + 2fi 4 Y 4 Yi 2 + 2fi b Y 5 YJ ) 2 )dm. 
2vr J S 2 2 2 

A direct computation shows that 

/ Y?dm=^, [ Y 2 2 dm=^, / Y 3 2 dm=^, 

U 1 15 y §2 2 45 y§2 3 15 

which implies that 

4 5 

( A + 6 -^)E^ = - 

b i=l 

Hence, A = ^ — 6. Specifically, Gh has a positive eigenvalue if and only if a > 7.5. □ 
When /i = h^ n^i = 1,2), the problem becomes more complicated. Kuzzu and Doi [9] 
conjectured that all the eigenvalues of Gh are non-positive, and Ker Qh = {0 • IZh, 8 6l 3 }. 
Here we will give a rigorous proof of Kuzzu and Doi's conjecture when h is a stable critical 
point. Let us introduce an important operator Ah defined by 

A h cj) = f -K ■ (hK<t>) 
The operator Ah has the following properties: 

Lemma 4.1. The operator Ah is a one-one mapping from H$(§ 2 ) to V (S 2 ). We denote by 
A^ 1 its inverse. Then it holds that 

Ah = A* h , {A h <t>A)>^ (A^<f>,<f>)>0. 
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Now we introduce another important operator T~Lh denned by 

n h f A =t+uf. 

We have the following important relation: 

Q h f = -AhHhf. (4.4) 
Basically, we can reduce the spectral analysis of Qh to that of T-Lh- 

Proposition 4.2. If h is a critical point of A[f], then QhAh is a symmetric operator and 

Moreover, if h is a stable critical point of A[f], then QhAh is a non-positive operator, and Qh 
has only non-positive eigenvalues. 

Proof. The identity follows from (|4.4p . Since h is a critical point, we have by (|3.ip that 

A h <t> = -Ml ■ K(j) - Kh ■ K<j) = -hK ■ K4> + hK{Uh) ■ He/). 
Then for any ip, (p £ H 2 (S 2 ), we have 

(g h A h ip, 4>) = -{KA h ijj + hlUAAh^) + A h ipKUh,K(t)) 

= -(KA h i/) + A h i>TUAh,Tl<t>) + (UA h ip,TZ- {hn<j))) 
= (A h ip, K-K4>- K(Uh) ■ K4>) - (UAhtp, A h <f>) 

= -(A h ^,^-)-{UA h i>,A h <j>). 



Specifically, 



(G h A h <j>,<p) = -(A h 4>, ^ + UA h <t>). 



h 

This means that QhAh is a non-positive operator if h is a stable critical point. Furthermore, 
if 4> is an eigenfunction of Qh associated with the eigenvalue A, then we have 

> {QhAhA^^A^cf) = (Qh^A^cP) = X^A^cf). 

Hence, A < 0. □ 
Now we establish a lower bound of the operator T~Lh- 

Proposition 4.3. Let hi = h„ ljn . For any f 6 Vq, there holds 

(H h J,f)>0, 

and the equality holds if and only if f € span{7£j/ii, i = 1,2,3}. Moreover, there exists a 
positive constant cq depending only on r\\ such that if f satisfies 



[ /(m)i, 1 K/ ll dm=0, 

J§ 2 



then we have a lower bound 

(H h J,f)>co(fJ). 
We need the following key lemma. 
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Lemma 4.2. Let r/ = rj(a) be determined by 113.2$ . and define A^irj) = z k e r)z2 dz. Then 
there hold 

A k+2 = — -(k + l)^, A = a(A 2 -A 4 ). 
rj lr\ 

Moreover, Sc g£ > when rj > rf ; d °Q^ < when r/ < rj* . 

Proof. The first equality can be easily verified by integrating by parts. While, the relation 
(|3.2|) is equivalent to 



6ae v - (3 + 2ri)aA = Ari 2 A <^=> A = a(A 2 - A 4 ). 

dr/ 



In order to prove the second statement, it suffices to show that the equation d"^) = has 



only one root, since - = 0. We have 




— (e"" {A (A 4 - A 6 ) - A 2 (A 2 - A 4 

= 1® f 1 f\ x ^ + x i y 2_ x 6_ y 6 + x i + y^_ 2x 2 y 2 )eV ^W~l) dxd y 

2 on J _i J _i 

( x 2 _ y 2 )2(1 _ x 2_ ^^(zW-l)^ < 0. 



7 1 

1 J-l 



Hence, Aq{A 4 — Aq) — A 2 (A 2 — A 4 ) = has only one root. Then from the fact that 
da(n) _ / A y _ A 2 (A 2 -A 4 ) - A (A 4 - A 6 ) 



dr, \A 2 -A 4 J {A 2 -A 4 y 

we know that — = has only one root. □ 

Proof of Proposition 14.31 Without loss of generality, we may assume n = (0, 0, 1). Intro- 
duce the sphere coordinates (#,</>) £ [0, 7r] x [0, 2tt] with m = (sin cos 0, sin # sin 0, cos#). 
Hence, 

g r;(cos 6) 2 

We make a Fourier expansion for / with respect to the variable </>: 
f = a (6) + Y, M#) cos (^) + W) sin(A^)) . 

k>l 

Noting that the area element dm = sin 9d9d(j), we make a change of variable z = cos 6 to get 
= 2^ ( f f'Az) ■ ( /' e-* 2 {24 + £>I + b 2 t )}dz) . 



fc>l 
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Routine computations show that 



(Uf,f)=-af [ (m • m') 2 /H/(m')dmdm' 

3 . 

a ( / Triimj f (m)dta 



= -2aTr 2 {(J (1 - z 2 )a dz) 2 + 2( J z 2 a dz) 2 + ( J zy/l - z 2 ai dz) 2 
+ ( J' zVl-z 2 hdz) 2 + -^f 1 (1 - z 2 )a 2 dz) 2 + ^(J 1 (1 - z 2 )b 2 dz) 2 y 

We use Lemma 14.21 and Cauchy-Schwartz inequality to get 

2vr 2 ( J 1 e^ds) • ( ^ e-^ 2 a 2 dz) -2c^ 2 ( ^ z^l - ^ 2 «idz) 2 

= 2a7r 2 (^ 1 e r ' z2 z 2 (l - z 2 )dz) • ( ^ e"^ 2 a 2 dz) - 2avr 2 ( ^ z\J 1 - z 2 a x dz^ > 0. 

Moreover, the equality holds if and only if a\{z) = Ce vz2 zy/l — z 2 or for some constant C. 
If ai(z) satisfies f_ l e vz zy/l — z 2 a\(z)dz = 0, then 



jf e"*V(l - z 2 )dz) • ( J e"^ 2 a 2 dz) - ( J zVl-z 2 ai d z y 
= (jf e"*V(l - z 2 )dz) ■ ( J e-^ajdz^j - ( J z y/l - z 2 (l - Ceroid. 

> ( jT 1 e ^ 2 (z 2 (l - z 2 ) - z 2 (l - z 2 )(l - Ce^) 2 )dz) • ( j' e^ajdz) 
= (^ 1 e 2 ^ 2 z 2 (l-z 2 )(2C-C 2 e^ 2 )dz) • ( j' e'^ajdz) 

> c ( V ) J afdz, (4.5) 



if we take £ small enough. We denote 

W^rj) = J e^ 2 (l - z 2 )(hz 2 - l)dz, 
which is positive for r/ > by noting that 



f 1 

y (e^ 2 -e^)(l -z 2 )(5z 2 - l)dz > 0. 
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By Lemma 14.21 and Cauchy- Schwartz inequality, we get 

2^ 2 ( J' e**dz) • ( ^ e^ z \ldz) - ^ 2 ( J* (1 - z 2 )a 2 dz 

= 2aTr 2 (j\ 1 > z2 z\l - z 2 )dz) ■ ( J' e-^ 2 a|cb) - ^a^ 2 ( j' (1 - z 2 )a 2 d; 

= l a7! 2^j\^^ _ z 2 } 2 d ^ . ^ J' e -^ 2 2 d ^ _ l ayr 2^ ^ (1 _ z 2 )(l2dz 

+laTT 2 Wi(ri) J e'^aldz^coir)) ^ a\dz 

for some c^irf) > 0. 

In the following, we take rj = f]i(a) for a > a* . By Lemma 14.21 we have 

da(y) = A 2 {A 2 - A A ) - A (A 4 - As) = 3A 2 + 2A A 2 - 5A A 4 
di] (A 2 -Atf 2r] (A 2 -A 4 y 

which implies that 

3(A A± - A 2 ) < 2A (A 2 - At). 
Then using the fact a$dz = and Cauchy-Schwartz inequality, we infer that 



J (1 - z 2 )a dzy + 2^J z 2 a dz 

= 3 (^-i)/>^ 

< 2{A 2 - At) J e^cgd* = ^(J ^dz 



e ^ 2 agdzp, 



which implies that 



2 

a 



( y e^dz) • ( y e-^agds) - (J (1 - z 2 )a dz) - 2^j z 2 a dz) 

>Co(r?) y Ogdz for some cq{t]) > 0. 
Summing up all the above estimates, we conclude that if rj = rj\(a), then 

(n hl f, f) = (j^—, f) + (uf, f) > coin) £ f 1 a 2 dz, 



and hence, 

/ / 



f) + (Uf,f) = 
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if and only if / is of the form C\Z\f \ — z 2 e nz cos (ft + C^-zVl — z 2 e r,z sin</>, which belongs to 
spanjT^/i^n, i = 1,2,3}, since we have 

Kh run =2? ? (m x n)(m-n)e , ' (m n)2 = 2r/( sin (9 sin 0, - sin (9 cos 0, 0) cos 0e v cos2 9 

=2ry( sin^, — cos (ft, 0) z\J 1 — z 2 e vz . 

This proves the first statement of Proposition 14.31 To obtain a lower bound of , we 
decompose g into 

f = fi + h, fi € span{^/ii,i = 1,2,3}. 
If <? satisfy J" §2 f(m)-A^TZhidm = 0, then we have 

/ / 1 (m)^ i 1 / 1 dm = - / / 3 (m)^>dm, 

which implies that 

(/i,/i)<C(/ 2 ,/ 2 ), 
since ~ (fx,A^ fi). This together with (|4.5p and (|4.6p gives 

(7^/,/) = (n hl h,h) > c (/ 2 ,/ 2 ) > c (/,/). 

The proof is finished. □ 
We define KerS S n t f e ij2(§2) . = }. 

Theorem 4.1. Let hi = h Vi;n ,i = 1,2. For a > a*, it holds that 

1- Qhi has no positive eigenvalues, while Qh 2 has at least one positive eigenvalue; 

2. KerCy/jj = |0 ■ !Zh\; G R 3 } is a too dimensional space; 

3. If (ft £ KerQ^, then T-Lhj(ft = 0; 

Proof. Since /ii is a stable critical point of A[f], Gh 1 has no positive eigenvalues by Proposition 
14.21 From the proof of Proposition 14.31 we know that there exists g G Vo(S 2 ) such that 

(n h2 g,g)<0. (4.7) 

Assume that all eigenvalues {Xk} of Qh 2 are non-positive. We denote by Ek the eigen- 
subspaces of Qh 2 corresponding to Xk- Then for iftk G Ek,ifti G Ei(k ^ £), we have 

Afc^fc,^ 1 ^} = {GiftkiA^t) = (Gifti,A^ift k ) = Xi(ift k ,A^ifte)- 

Hence, (iftk, A^ifte) = for k ^ £. We write g = J2 k tftk with iftk G Ek- Then 

(n h2 g,g) = -(Gg,A£g) = - X k (ift k ,A^ift k ) > 0, 

k 

which leads to a contradiction with ()4.7p . Thus, <7/i 2 has at least one positive eigenvalue. 

If (ft G Ker^/jj, then = constant. Hence, 0) = 0, and then G spanj^/i!, i = 

1, 2, 3} by Proposition 14.31 On the other hand, if (ft = TZhi, then we find 

n hl (ft = 7£(ln/ii +Uhx) = 0. 

This proves Ker^/ ll { ■ !Zh\; G R 3 } and the third point. Due to n • !Zh\ = 0, KeiQ^ 
is a two dimensional space. □ 
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Finally let us give a characterization of the functions in KerQ^ n , see also [9] . 

Proposition 4.4. If ifto G KerC/^ n , then ifto takes the form (9, (ft) 

V>o = © • e^otcostf), 
in the spherical coordinate, where = (— sin^>, cos (ft, 0) and go satisfies 

1 d / . „d5( \ c/o du d5o du , . 



sin0d0V d0/ sin 2 d0 dfl " d9 



Proof. Note that KerC/f = A h 1 KerC//, . Hence, ibn E KerC/f if and only if there 
exits a vector such that 

which is equivalent to 

K-Kipo -Ku Q -Kift = -0 -Ku , (4.9) 

where uq = Uh Vl)T1 is a function of m • n. We take 6 be the angle between m and n, and 
rewrite (|4.9p in terms of in the spherical coordinate (9, (ft) as 

1 8 ( . Q dift \ 1 <9 2 V>o du dift du 



sm9 89\ 89 J sin 2 9 deft 2 d9 89 9 d9 ' 

We rewrite ifto(9, (ft) as 

^{9, (ft) = • e^go. 

Then it easy to find that go satisfies (|4.8p . □ 

5. Lower bound of a bilinear form for the linearized operator 

In the inhomogeneous case, the linearized operator Q e h around h is given by 

Q e h f = Tl- (TZf + KRUef + fRUh) . 

To justify the small Deborah number limit for the inhomogeneous system, the main difficulty 
is that the elastic stress in the velocity equation is strongly singular(a loss of -). To overcome 
it, we have to establish a precise lower bound for the following bilinear form: 

(G £ h f,uy), nf= f - + u £ f. 

When £ / 0, the orthogonal structure is destroyed such that the interactions between the 
part inside the kernel and the part outside the kernel of Qh become very complicated. We 
find a coordinate transformation and introduce a generalized kernel space of such that the 
interactions between two parts can been seen explicitly, then a lower bound is obtained by 
very subtle calculations. 

5.1. New coordinates frame. At each point x, we choose a right hand cartesian coordi- 
nate frame (ki(x), k2(x), ka(x)) such that k3(x) = n(x), and ki(x), k2(x) depend on n(x) 
smoothly. For instance, under the assumption that |ni(x)| < 1 — cq for all x 6 O, we can 
take 

n x (1,0, 0) T ( n 3 n 2 n t 

k 2( X ) = 1 n nNTl = > 



nx(l,0,0) T | V ' ( n 2 + n 2)l/2' ( n 2 +n 2)l/2 

, ^ /, 2 , 2^1/2 n i™2 rnn 3 

k 2 x n = ((n 2 + «,)', ~ (re 2 + n§)1/2 , - (n 2 + n§)1/2 
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At each point x, let (9, 0) be the sphere coordinate on the unit sphere § 2 , that is, 

m = sin 9 cos 0\ii (x) + sin 9 sin 0]s.2 (x) + cos 9h-^ (x) 
= A(x) • (sin 9 cos 0, sin 9 sin 0, cos 9) T , 
where the matrix A = [ki k2 k3]. We set 

= — sin 0\l\ + cos 0^2-, = —(cos 9 cos <£ki + cos 9 sin <£>k 2 — sin 6\l^). 

We denote rh = (sin0cos(^,sin^sin(^,cos^) r , hence m = A(x) • rh. 

In this coordinate, the rotational gradient operator 1Z = m x V m can be written as 

9 - 1 9 

1Z = ( — sin <£ki + cos 0k.2 ) — - — (cos 9 cos 0ki + cos 9 sin (^k2 — sin (^3) 



39 sin <9£ 

. / « d „cos9 9 ^ d „cos9 9 9 n t , N 

= A • - sin 09 — 7 - cosw -— — ■ ,cosm — - — sm<2 -TrT,Trr • 5.1 

v 89 sinO d<p 86 sin 9 dip dip 1 

We also have 

_ „ _ / df 1 df-. , dg 1 dg s 

TZj ■ Tig = (e^— + e„- - — ) • (e^— + - — ) 

y V d9 e sin9dip J v v d9 d sin9 9ip ! 

= 9f9g^ I df dg 
99 99 sin 2 9 90 80' 

5.2. The Maier-Saupe space and the lower bound inequality. For any / 6 L 2 (fi x § 2 ) 
with Jg 2 /(x, m)dm = 0, we decompose it as 

/ = ao( x > 8) + («fc( x 5 cos fc<£ + 6fc(x, 9) sin fc<£) , (5.3) 
fe>i 

with Ofc(x, 0) = afc(x, 7r) = 6fe(x, 0) = 6fc(x, tt) = for A; > 1. We set 



A k = e r > cos29 (cos9) k sin9d9, f o 0) 
Jo 



^rj cos 2 I 



2nA Q 

We further decompose the coefficients ao,a\,a2, &i, 62 as follows 

a (x, (9) = Co(x)/o(^)(cos 2 (9 - A 2 /Ao) + 7o(x, 0), 
ai(x,(9) = Ca,i(x)/o(0)sin0cos0 + 7ai i(x, 0), 
61 (x, 0) = C M (x)/ ((9) sin cos (9 + 7M (x, 0), 
a 2 (x, (9) = Ca, 2 (x)/ o (0) sin 2 + 7a , 2 (x, (9), 
& 2 (x,0) = C6,2(x)/oWsin 2 ^ + 76i 2(x,^), 
where the functions 7o,7a,i> * * * j7j>,2 satisfy 

/ 7 O (x,0)sm0d0 = 0, / (3 cos 2 9 - l)7o(x, 9) sin (9d(9 = 0, 
io io 

/•7T rn 

/ sin0 cos 7a) i(x, 9) sin0d0 = 0, / sin 6* cos 07^1 (x, 9) sin0d0 = 0, 
Jo J 

/*7T /*7T 

/ sin 2 07 a>2 (x,0)sin0d0 = 0, / sin 2 07& >2 (x, (9) sin0d0 = 0. 
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Noting that 

I /o(cos 2 0- J 4 2 Mo)(3cos 2 0-l)sin^ = ^-(A4-4^) > 0, 
Jo ^0 ^0 

hence Co is uniquely determined. Obviously, Ca,i> • • • , Cfe,2 are also uniquely determined. Thus, 

the above decompositions make sense. 

The space spanned by the following five functions 

/o(#)(cos 2 — A 2 /A ), /o(#) sin# cos # cos /"oC^ 1 ) sin ^ cos 6 1 sin <^>, 

/o(#)sin 2 (9cos(2<£), and f (0) sin 2 9 sin(2^) 

can be viewed as a generalized kernel of the operator Qf , and will be called as the Maier- 
Saupe space. The kernel of Qh is spanned by the second and the third function. 
We denote 

M(x) = / mm/dm, M(x) = / mm/dm = / (A • rh)(A • rh)/dm, 

Js 2 Js 2 Js 2 

Nij (x) = A ki A tj (g £ * M kl ) , iV (x) = 2iy 33 - N n - N 22 , iV 2 (x) = N u - N 22 . 

Proposition 5.1. (Lower bound inequality) Let h n M jTI be a stable critical point of A[f]. 
Then there exits c > such that any f £ H 1 ^ x § 2 ) with L 2 /(x, m)dm = 0, there holds 

(SLJ> J) >c{ f r (To' + 7a,l + ill + 7a,2 + 76,2 + J>* + S ' 1U 

+ [ ((2Co - aiV ) 2 + (Ca,i - 2aiVi 3 ) 2 + {Cb,i ~ 2aN 23 ) 2 + (Ca, 2 - aN 2 f + (4,2 - aN 12 f 
Jn 

Remark 5.1. This inequality gives a good bound for the part outside the Maier-Saupe space, 
and a weak bound for the part inside the Maier-Saupe space. 

Remark 5.2. By letting e tend to zero, the lower bound inequality implies that 

(&»,„/, W) > c(n fj), n = n haiV , 

which can also be deduced the following simple argument. Noting that T-L(KeiQh av ) = 0, we 
may assume that f £ (KerQ^ ^)- L . Then by Poincare inequality and Proposition ^.!^ we have 

(Hof,AHaf) > (KHof,KHof) > (H f -Wf \H f -Wf) 

v2 



> ( ^° / (/ ^ /,// > c{H f, f) > c(f, /), 



where H f = ^ J§ 2 "Ho/dm. 



5.3. Proof of the lower bound inequality. First of all, we can get by direct calculations 
that 

U £ f = a / g £ (x — x')(l — mm : m'm')/(x / , m'jdm'dx' 
Jn Js 2 

= —a mm 



: / 5£ (x-x')M(x')dx' 
Jn 

a(A-m)(A-m): / g £ (x - x')M(x')dx' 
Jn 



-omiimj 



AkiAij(g £ * (Aki'Aij'Mi'f)). 
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For the simplicity, we denote fo = h Tun in what follows. Using (|5.ip - (|5.3p . we get by very 
tedious calculations of competing the square that 

(QhJ^fJ) = </o^(y- +UJ),K(l- +UJ)) 
Jo fo 

f f 
= (foK(— - arhirhjNij),lZ(— - arhirhjNij)) 

Jo Jo 

= f f fo\do{^- - Qm i m J iVj J )| 2 dmdx+ / / ^° - |<9y(^- - arhirhjNij) | 2 dmdx 
Jn Js 2 Jo Jn J§ 2 sur 8 Jo 

= 11 /o|^(?) + «(2iV 3 3-iVii-iV 22 )sin0 cos 0| 2 drhdx 
Jn Js 2 JO 

+ 111 Md§&) -2aN 13 cos 28 \ 2 + J^(^±-2aN 13 sm§ cos 0) 2 dmdx 
2 Jn Js 2 Jo sin z 8 Jo 

+ 111 fo\dA-2aN 23 cos28\ 2 + ^^(^-2aN 23 sm8cos6) 2 dnidx 
2 Jn Js 2 Jo sin^ 8 Jo 

+ \ I I /o|^(^)-a(^ii-^ 22 )sin^cos0| 2 + -^(^-a(^ii-^22)sin 2 ^) 2 dmdx 
2 Jn Js 2 Jo sin z 8 Jo 

+ 111 fo\ d - 2aAf i2 sm 8 cos 8? + -%(^- aN 12 sin 2 tffdmdx 
2 Jn Js 2 Jo sur 8 Jo 



2 ^ A} 7§ 2 /o /o sin 2 8 fo 

Making a change of variable z = cos 8, we obtain 

f (l - z 2 )(d z (^) - a(2N 33 - Nu ~ N 22 )z) 2 dzdx 
nJ-i Jo 



+ 2 



\ [ ^ fo\d z (^)VT^+2aN 13 (2z 2 -l)\ 2 + r^{^-2aNi 3 zVT^) 2 dzdx 
2 Jn J -i jo 1 — z Jo 

+ 1 I r fo\dAVT^ + 2aN 23 (2z 2 -l)\ 2 + -^{^-2aN 23 zVT^) 2 dzdx 
2JnJ~i Jo l-z z Vo 

+ \ I I /o(l-^)(^(^) + a(iVii-iV 22 )^) 2 + -^(^-a(iVii-iV 22 )(l-z 2 )) 2 dzdx 
^ Jn J -i Jo 1 — z jo 

+ \t I'' Ml-z 2 )(d z ( h f) + 2aNi 2 z) 2 + -^& -aNi 2 (l-z 2 )) 2 dzdx 
* Jn J -i Jo 1 — z jo 

+?E/„£/"i 8 s(x)i 2 +/»i s »(x)i 2 + -4^^. 



2 7f7 7§ 2 /o fo sin 2 /o 

• Lower bound for the term including a 

To deal with the cross term, we need to introduce a slightly different decomposition 

a (x,z) = Co(x)/ (^)(^ 2 - 4^) +7o( x ^)> 
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where 70 satisfies 70 (x, z)(3z 2 — 1 — 2r\z 2 (\ — z 2 ))dz = 0. From the fact that 

j f (A z 2 - A 2 )(3z 2 - 1 - 2 V z 2 (l - z 2 ))dz 
= 3{A A 4 - A 2 ) + 2 V (A 2 (A 2 - A 4 ) - A (A 4 - A 6 )) > 0, 
we know that £0 are uniquely determined. Then we have 

/o(l - z 2 )\3 z (^) - aiVo(x)z| 2 dzdx 
Jo 



n j-i 



■■J J 1 HI ~ z 2 )(d z (^j^)) 2 dzdx+ ( J (2Co - aiV ) 2 dx) ( j\l - z 2 )z 2 f dz 

-2/ (2,4oCo(x)-aiVo(x)) / d z (^^)f (z)z(l - z 2 )dzdx 
Jn J-i Jo 

J f fo(l ~ z 2 )(d z (^j^)) 2 dzdx + ( J (2Co - aiV ) 2 dx) ( j\l - z 2 )z 2 f dz) 

+ 2 I (2A)Co(x) -aiV (x)) f 7o (x, z) (l - 3z 2 + 2 V z 2 {l - z 2 ))dzdx 
Jn J-i 

J f /o(l - z 2 ){d z (^j^)fdzd^ + ( J (2Co - aiV ) 2 dx) ( j\l - z 2 )z 2 f dz 



Recall that we have decomposed ao(x,z) as 



A 2 

ao(x, z) = ( (x)f (z)(z 2 - — ) +7o(x, z), 
thus (Co - (o)fo(z)(z 2 - j^) = 7o - 7o, which implies that 



Co(x) - Co(x) = J 7o(3z 2 - l)dz/(3A A 4 - 3A 2 ) 



Thus we have 



. 1 / (l-z 2 )(^(^)) 2 dzdx= / / f (l- z 2 )(d z A-2A (Co-Co)z) 2 dzdx 
nJ-i Jo JnJ-i v Jo 7 

<c [ /' 1 / (i-, 2 )(a 2 (^)) 2 + 7o 2 ^dx. 

Jn J-i v Jo 7 

Then we infer that 

1 / (1 - z 2 )\d z (^) - aiVo(x)z| 2 dzdx 
Jo 

{ J f /o(l - z 2 )(d z (^j^)) 2 dzdx + J (2 Co - aiV ) 2 dx |\l - z 2 )z 2 f dz} 

>c{ f [ 7o 2 + (^ 7o ) 2 d2dx+ / (2Co(x)-aiVo) 2 dx), (5.4) 
l JnJ-i Jn } 

where we used the following Poincare type inequality in the last inequality: 



'fi J-i 

> c 
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Lemma 5.1. There exists c > such that if j(z) satisfies 



J jdz = 0, J (3z 2 - l)7<Le = 0, 



then there holds 



C f (l-z 2 )(d z (^)) 2 dz>c [\ 2 dz. 

J-l JO J-l 

Proof. We define 7(111) = 7(111 • n). By the assumption, we know 

/ mm7(m)dm = 0, ^07=7-- 
Js 2 Jo 

Hence, 7 G (Ker£?y )^ and we have 

/ /o|7e^ /o 7| 2 dm = 2vr C / (1 - z 2 ) (d 2 (-^)) 2 dz. 

J§ 2 J_x Jo 

Set C = f iHf jdm. It follows from Poincare inequality and Proposition I4.3I that 

[ tvou -|2j >> /" /oy - (/§2(^/o7-C')7dm) 2 
/ fo\KUf j\ dm > c / (H/ 7 - C) dm > c— 
is 2 J§ 2 



is 2 7 2 d 

. (/§ 2 (^/o7)7dm) 2 f _ 2 2 

> c „ — > c 7 dm = c47r / 7 dz, 

Jg 2 7^dm 7s2 



m 

1 



1 



which completes the proof. 

• Lower bound for the terms including a±,b\ 
Recall that we have a decomposition for a±(x, z) as 



ai(x,z) = C,a,i{-x)fo{z)z\/l - z 2 +7 a> i(x,z), 
where §_ x 2VI — 2; 2 7 aj idz = 0. Thus, we can get 

/ t ~ 2aN 13 zVT^) 2 dzdx 

Jn J-i 1 - z Jo 

> I T /o P' 1 ^ + (Ca,l(x) - 2aAT 13 (x))zv / r^ 2 ) 2 dzdx 
Jn J-i v Jo 7 

= f C 7a2 ' l( „ X ' Z) dxdx+ / (Ca.i(x)-2aiY 13 (x)) 2 dx [' f z 2 (l - z 2 )d 
Jn J-i Jo Jn J-i 

A lower bound for the terms including h\ can be obtained in the same way. 
• Lower bound for the terms including 02,62 
We have decomposed 0,2(3:, z) as 

a 2 (x, Z) = Ca,2(x)/ (^)(1 - Z 2 ) + 7a,2(x, z), 
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where 7a,2(x, z)(l — z 2 )dz = 0. Then we have 

1 J^(^_ aiV2 (l_ z2) ) 2 dzdx 
1 - z 2 Vo 



n J-i 



> f f 1 4/ ( 7a ' 2( „ X ' Z) + (Ca, 2 (x) - aiV 2 (x))(l - z 2 )) 2 dzdx 
JnJ-i v Jo 7 

= / / 1 ^4^dzdx+ / (Ca, 2 (x)-«iV 2 (x)) 2 dx / (l-z 2 ) 2 dz. (5.6) 
JuJ-i Jo Jq J-i 

We can obtain a similar bound for the terms including o 2 . 

Finally, the lower bound inequality follows from (|5.4I) . (|5.5H and (15 .6p . □ 

6. Small Deborah number limit for the homogeneous system 

This section is devoted to justifying the small Deborah number limit for the homogeneous 
system (|2.ip . For the simplicity of notations, throughout this section we denote 

ll/IU = ll/lltffc(S*)> \f\v = \\f\\Lp(§i), (f,g)= /(m)j(m)dm. 

6.1. Hilbert expansion. As in the fluid dynamic limit of the Boltzmann equation [T], we 
make the Hilbert expansion for f £ (m,t): 

3 

f(m,t) = Y,e k fk(m,t)+e 2 f e R (ui,t). (6.1) 

fc=0 

Plugging it into (I2.ip and collecting the terms with the same order with respect to e, we find 
that fk(m,t)(k = 0,1,2,3) satisfies 

Kfo + foKUh = 0, that is / (m,t) = h a[t) (m), (6.2) 
^ = 0n(t)fi ~ n ■ (m x k • m/o), (6.3) 
^ = a„ ( t)/2 + ft • (/«) - ft • (m x K • m/i), (6.4) 

= ^n(t)/3 + ft • (/iftW/ 2 ) + ft • (/ 2 ftW/i) - ft • (m x k • m/ 2 ). (6.5) 

The global in time existence of the Hilbert expansion is nontrivial, since f\ satisfies a 
nonlinear equation. However, we find that the part of f\ inside the kernel of £7 n (t) satisfies a 
linear equation by Lemma [ 



Proposition 6.1. Let n(t) be a solution of \2. 3\) on [0,T] with A given by A2.5\) . We can 

construct smooth functions /fc(m, t)(k = 0, 1, 2, 3) £ Vq defined on [0, T] suc/i £/ia£ i6.2\) - (675\) 
hold on [0,T]. 

We need the following two lemmas in order to prove it. 
Lemma 6.1. [HE] Let /o(m,i) = /i n ( t )(m). Twen /o(m, i) satisfies 

(^+ft-(mx K -m/ ),v)=0 (6.6) 
/or any ^ G Ker^*^ i/ and onh/ if n(t) is a solution of \2. 3\) with A = A(a). 
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Lemma 6.2. For any (ft, (ft € KerG n , there holds 

<ft = -KU<ft, ((Ucft) 2 ,0) = O. 

Proof. For (ft G KevQ n , there exists £ R 3 such that = 0- 1Zh n by Theorem 14. II Due to 
ln/i n = —Uh n , we see that 

7£/i n = —h n 7ZUh n = —h n U!Zh n . 

Hence, (ft = —h n U(ft. To prove the second equality, we choose the sphere coordinates such 
that n = (0, 0, 1) and m = (sin 9 cos tp, sin 9 sin tp, cos 9). Let (ft = • 7£/i n for some S R 3 . 
Then we have 

eft = r/Zin cos 9 sin # (Gi sin tp — @2 cos <£>) , 
^ = V^n cos sin 9 (0i sin 99 — ©2 cos tp) . 

It is easy to check that 

I (0i sin <p — 02 cos (^) 2 (Oi sin 93 — 02 cos <p)d(p = 0, 
Jo 

which implies that J g2 p-^dm = 0, or equivalently {{U(ft) 2 ,(ft) = by = —h n U(ft. □ 

Proof of Proposition 16.11 Let us first solve f\ and write f\ = (ft{t) + it), where (ft £ 
Ker^ n , € (Kert/*)- 1 -. Then will be determined by (|6.3p . while will be determined by 
(|6.4p . However, (|6.3j) has a solution if and only if 

(?j£+K-(mXK-mf ),il>) =0 

for any ^ £ Ker<7*. From Lemma l6.lt this is equivalent to require that n(t) is a solution of 
(|2.3|) with A = A (a). Given (ft- 1 -, (|6.4p implies that satisfies 

(?±^) = (K-(((ft + ^)KU{4> + <ft X )) - Tl ■ (m x (k ■ m)(0 + <^)) , V> - (^f^) 

for any ^ € Ker<5*. Since KerC?* is a two dimensional space, we take ifti,ip2 as a base of 
Ker<5* and write (ft = a\{t)ifti + 02(^)^2- Then we can get an ODE system for (a\(t), 02 (i)) ■ 
For any tp € KerC?*, we write ip = A^ 1 ^ with (ft 6 KerC/. Due to Lemma 16.21 we find that 

{n-{(ftnucp),i}) = (tz ■ {(ftWA^),A^(ft) = -{(jmucft^A' 1 ^) 

= {Ku^nu^izA^ 1 ^) = i<^[(^) 2 ]> n ^ n 1 0> = ^((UcftfA) = 0. 

This means that (ai{t), 02 (^)) satisfies a linear ODE system, hence is global in time. 

Once fi is determined, we can get /2 and by solving (|6.4p and (|6.5p in a similar way 
(note that the equation for fi is linear). □ 

6.2. Error estimates. This subsection is devoted to proving Theorem l2.11 Due to the weak 
nonlinearity of the kinetic equation (|2.ip . given the initial data /g € i/ 1 (S 2 ), it is easy to 
show by standard energy method that there exists a unique global solution f £ (m,t) to (12. ip 
such that 

f £ e C([0,+oo); J ff 1 (S 2 ))nL 2 (0,T; J ff 2 (§ 2 )) for any T < +00. 
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Thanks to Proposition 16. II and (|2.ip . it is easy to find that f R satisfies 

J^(m, t) =- £ G n f £ R -K-{ m xK. m f R )+eK- {f R WAf R ) 

3 

1=1 

where A is given by 

l<*J<3,i+j>4 

To complete the proof, it suffices to prove that 

||/r(*)II-i ^ C fora11 0<t<T. 
For this purpose, we need the following lemmas. 
Lemma 6.3. 



■at' n j^— - n v at 

Proof. Direct calculations give that 

_ _ | R . { kk s ) = -n • - n . (a.r| 

Replacing 5 by .A^ 1 *?, the lemma follows. 

Lemma 6.4. For any vector field V G C 1 (S 2 ), i/iere /10/ds 

(ft- 0*7), ^n 1 /) < cflyioo + IwiooX/,^- 1 /). 

Proof. Let V = {V u V 2 , V 3 ) T ', ft/ = (i?i/, i? 2 /, i? 3 /) T , and 5 = ^/- Recalling A 
-ft ■ (h n Kf), we get 

(ft ■ (V f),A~ x f) = - (RilViR^R^g) 

= - (h n R j9 , (RjV^Rig + ViR 3 R t g) 

= - (KR 3 g, {RjV^Rig - Vie kji R k g) + (h a R j9 , V^R.g) 

= - (h n R j9 , (RjV^Rig - Vie k ^R k g) - ^(R t (h n Vi), (R^g) 2 ), 
which implies the lemma by using the fact that 

\n g \ 2 <c\\f\u<c{f,A- 1 f). 

The proof is finished. 
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Now we are in position to prove (|6.8p . With the help of Lemma 16.31 we make the energy 
estimate for (16.71) to obtain 



1 



n JH/ 2\ n V ^ 



i=i 



Since /i n is a stable critical point, we know from Proposition 14.31 that 

(Gnf R, A' 1 f R ) = -(H a f R ,f 

We infer from Lemma 16.41 that 

-(K- (m XK-mf^^^fn) < C{f R ,A n l f R ), 

{eK-{f R 1ZUf R ),A- l f R ) < Ce{f R ,A- l f R )K 

<^-(/^/ 4 ),^ n 1 /^><C||/,||_ 1 </^^ n 1 ^>, 
while the other terms on the right hand side are bounded by 

3 1 

(E eJ " 1 |^l2 + |9 i /o| 2 )</^-An 1 ^>+ePI|-l</|j,-4n 1 /^ 5 . 
i=l 

Here we use the fact that \1Z k Uf R \oo for k € N can be bounded by (f^A^ 1 /^) 2 and 
i 

(ffrAn 1 /^) 2 ~ ||/^||_i. In conclusion, we obtain 

^\\ffl\±i<c{\\mli + 4ffl\li + 4ffl\-i)- 

This implies (j6.8f) . □ 

6.3. The Lesile stress and coefficients. This subsection is devoted to proving Theorem 
12.21 We introduce the 2-order tensor Q2U] and 4-order tensor Q4[/] as follows 

Q 2 [/] = (mm-il> / , 

Q4a/37/x[/] = (rnampm^m^ - ^(m a mp5 lfl + m 7 m^5 a ^ + m a m 7 5 /3M + mpm^5 a7 
+m a m fl 8i3 1 + mpm^Sa^) + -^(8^8^ + 8^8^ + 8 a/J ,8^)^ . 

Lemma 6.5. It holds that 

Q 2 [/ i „] = (P 2 (m.n)) fen (nn-l), 

8 cry 8 p fi 
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Proof. We only prove the first equality, the proof of the second equality is similar but 
more complicated. Since both sides of the first equality are tensors, they are coordinate- 
independent. So, we may choose the sphere coordinates such that n = (0, 0, 1) and m = 
(sin 9 cos </?, sin 9 sin tp, cos 9) . Since h n (m) depends only on cos#, it is easy to check that 

/ mimjh n (m)dm = for i ^ j, 

Jn 2 

(m| - i)Mm)dm = ^(3cos 2 # - l)/t n (m)dm = ^<P 2 (cos 6)) h ^ 

( ( m i - ^n(m)dm = / (ml - -)h a (m)dm = -^(P 2 (cos 6>)}, , 
Js 2 3 J S 2 3 3 

which give the first equality. □ 

Lemma 6.6. Let f £ (t) be given by Theorem \2.1\ and P(x) be a smooth function on R. Then 
we have 

\(P(m.n(t))) fe[t) -(P(m.n(t))) hn(t) \<Ce. 
Proof. By the definition, we get 

<P(m • n)) /e - (P(m • n))^ = / P(m • n) ( £ e k f k (m, t) + e 2 f R (m, f))dm. 

^ 2 k=l 

Using the facts that /i, /2, /3 are bounded and 

/ P(m • n)/|(m,i)dm = / A n P(m ■ f £ R (m,t)dm 
Jn 2 J§ 2 

< p n P(m-n)|| ||/^||_i, 
the lemma follows. □ 
Now we are in position to prove Theorem 12.21 Direct computation shows that 

^Q 2 [/ £ ] =- / (mm - • ( - e ^(m x K ■ m/ c ))dm 

at e J S 2 3 

= -(m x 7£/i e m + mm x IZfi 6 ) fe — (2D : (mmmm) fe 

— D ■ (mm) fe + ft • (mm) fe — (mm) j= • (D + ft)) . 
So by Lemma 16.61 the stress a £ can be written as 

a 6 =— D : (mmmm)j£ (mm x TZfi[f £ ])fe 

2 s 

=— D : (mmmm)p — — ^2D : (mmmm) fs — D • (mm)^E 

+ ft • (mm) f e - (mm) f e • (D + ft) + ^Q 2 [/ £ ]) 
= ^D : (mmmm) ftn - t^2D : (mmmm), )n - D • (mm) ftn 

+ ft • (mm) hn - (mm) hn • (D + ft) + ^Qifin 



^[^ efc | /fc(m ' t)+e2 | / ^ (m ' t) 
k=l 



+ Ce. 
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For any constant vector U, V, we have 

e 2 U T • Q 2 [^/£(m, t)]-V = e 2 J (m -U)(m- V) [e^g n f R -%x«. mf R ) 

3 

• (f R KUf R ) + ^ e*- 1 ^ • (/f TWf R + KW/f ) + eA] dm, 
i=i 

which is bounded by C£||/ft||_i by using the argument of integration by parts. Then we infer 
that 

a 6 =^D : (mmmm)^ - - ^2D : (irimmm)/ )n - D • (mm)/j n 

+ ft • (mm) ftn - (mm) ftn • (D + fi) + ^ChM) + Ce. 
We see from the definition of Q4 that 

D : (mmmm)^ = D : Q 4 + : (mm) hn I + -(D : Q 2 + Q 2 : D) + — D. 

7 7 15 

Using Lemma 16.51 we can calculate 

a e =-(- 5 4 (D : nn)nn-y(D : nn)I - 5 2 (nN + Nn) 

8 10 2 ^2 \ 

+ ( 15 " 21 52 " 35^ D + + 7^ ( nn • D + D • nn)J + Ce. 



So, we conclude that 



cj £ - ^ - pi = ( n(-N - D n) - (-N - D • n)n + Ce = Ce, 

here we used the equation (12.3|) . This completes the proof of Theorem 12.21 □ 

7. Small Deborah number limit for the inhomogeneous system 

This section is devoted to justifying the small Deborah number limit for the inhomogeneous 
system (|2.8p - (|2.9p . Throughout this section, we will use the following notations. (,) denotes 
the inner product on L 2 (n x S 2 ) or L 2 (0). We also denote \\f\\ LP = ||/||LP(nxS 2 )(ll/llLP(fi)) 
for / defined on x § 2 (for / defined on 0), and || • \\jjo,k = |||| • l|.H*(n)||£2(s a )' 

7.1. Hilbert expansion. Due to the choice of g, we have the formal expansion to the 
operator U £ : 

U e f-Uf=f [ a|mxmf 9e (x-x')(/(x',m',t)-/(x,m',i))dm'dx' 
= / / a\m x m'| 2 5 (y)(/(x + y/ey,m',t) - /(x, m',i))dm'dy 

= [ f a\m x m'lVy) ( V ^-(y ■ V) fc /(x, m', t))dm'dy 
= [ [ a|mxm / |Vy)(V-^-(yV) 2fc /(x,m / ,t))dm'dy. 

We denote 

C4[/](x,m,t) = -l- / / a|mxmf S (y)(yV) 2f! /(x,m' 1 t)dm'dy for fc > 1. 
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Formally, we have 

Hef = U [f}+eU 1 [f}+e 2 U 2 [f} + --- , U [f]=Uf. (7.1) 
Then we make a formal expansion for the solution of (12.8p - (12.9|) : 



/ £ (x, m, t) = m, t) + e 3 f R (x, m, t), 

k=0 

2 

x,i) = ^2e k v k {x,t) +e 3 v £ R {*,t). 



V 

k=0 



Plugging them into f|2.8[) - (|2.9|) and collecting the terms with the same order with respect to 
e, we find that 

ft/o + foWJh = 0, that is, / = ^, n (t,x) ( m ); (7-2) 
and for the terms of order O(l), there hold 

^ + v • V/ =g fo h + ft • (/oW/o) - K • (m x (Vv ) T • m/ ), (7.3) 
dvQ 7 1 — 7 

— + v • Vv = - Vp + ^ Av o + ~2Jte V ' 1 ^ mmmm ) /o) 

-^{v-(mmx(ft/ 1 + /^A)>i + (V^i/o)/o}; (7.4) 

i+j + fc=l 

and for the terms of order O(e), there hold 

^-+v -V/i = S /o / 2 + ft- ( ^ fiKUjfk)-vi-Vfo 

i+j+k=2,j>l 

- ft • (m x (Vv ) T • m/i + m x (Vvi) T • m/ ) , (7.5) 
<9vi 7 1 — 7 / v \ 

-7^- + V • Vvi = -Vpi + ^ Av l + ~^T V ' ( 2^ D * : ( mmmm )/j J - V l • Vv 

-^{v-(mmx (ft/ 2 + £ f i KU j f k )) 1 + { ]T /iVE/,-/fc>i}; (7.6) 

and for the terms of 0(e 2 ), there hold 

^+v -V/ 2 = S /o / 3 + ft-( ^ /i^i/* 

i+j+fc=3,j>l 



- ft • ( ^ m x ((V Vi f • m)/,) - Vl • V/i - v 2 • V/ , (7.7) 



i+j=2 



+ v • Vv 2 = -Vp 2 + ttAv 2 + tjt^V • f 23 D< : (i 



dt * Re 2Re 



1-7 



{V • (mm x (ft/ 3 + ]T f i KU j f k )) 1 + < £ /iVE7,-/ fc }i} 

- vi • Vvi — V2 • Vvq. (7.8) 
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Proposition 7.1. Let (v ,n) G C([0, T];H 20 (£l)) be a solution of l2~JM - l2~Ti\) with A given 
by \2. 5|) and Leslie coefficients defined by \2. b)) - [27i\ ) on [0, T]. Then there exist the functions 
fi G C([0,T];i/ 20 - 4i (fi x S 2 ))(i = 0,1,2,3) and v< G C([0,T\; H 20 ~ Al (n))(i = 0,1,2) suc/i 
that (f7JP-(TOP ZioWs on [0,T]. 

We need the following lemma. 

Lemma 7.1. [5] Let /o = fr^nfx,*) ( m ) ■ ?7&en /o(x, m, i) satisfies 

(| + vo'V/o + K'(mx (Vv ) T • m/ ) + ft • (/ Kf/i/o), ^) L2(§2) = 
/or any ?/> G KerC/^ i/ and only if n(x, i) is a solution of A2.10\) with v = vo . 

Proof of Proposition 17.11 We denote by Pj n the projection operator from Vo(S 2 ) to 
Ker<5/- , and denote by P ou t the projection operator from Vo(§> 2 ) to (Ker^^)^. Let P- m fi = 
4>i j IPout/i — ipi- Now tpi will be determined by (j7.3|) . whose existence is ensured by Lemma 
17.11 Once ip\ is determined, it can be proved that the equation (I7.4D is equivalent to (12. lip , 
see [5]. Now we solve (0i,vi). In what follows, we denote by L((p,v) the terms which only 
depend on (</>, v)(not their derivatives) linearly. Let 0i = n ■ ft/o- We have 

(^ + v O -V)0! = (^ + v -V)n ± -ft/o + n ± -ft((- + v -V)/o). 
This means that 

4 +v »' v )*) = p »'( n± - K «l 



P 



out 



((^ + v -V)^) =P out (n ± -7e((^ + Vo -V)/ )) ^L(^), (7.9) 



P in ((^ + v • V)<k) = {-g t + v • V)0 X - L(&). (7.10) 

We also have 



VmiGfofi + K-ttiKUfa)) =0. (7.11) 

For a matrix k, we denote 

K{k) = P in (ft- (mx K-m/o)), £(«) = P out (ft • (m x « • m/ )), 
Binfai) = Pin (ft- {fonUxfa)), Boutifa) = F out (K- (foKUKh)). 

Taking Pj n for the equation of /i, we get by (|7.10p and (|7.1ip that 

(§- t + v o • V)0i = L(0i) - /C((V Vl ) T ) + B m (^) 

+ P in (ft • (faKU-ip! + i)<RU$i + i^iRUi)!)) - vi • V/o 
- P in (ft ■ (m x (Vv ) T • m/i)) + P in (ft • {foKU^)) + P in (ft • (/iftt/i/o)) 
= -/C((V Vl ) T ) + e in (^) + L(0i, vi). (7.12) 
Taking P out for the equation of /i, we get by (|7.9p that 
- (0/ o / a + ft • (faKUfa)) = + £((V Vl ) T ) - 

+ P out (ft • {fcTW-tyx + t/>iftZ#i + ipiRU^x)) 

- P out (ft • (m x (Vv f • m/i)) + P out (ft • (foKU^i)) + P out (ft ■ tfiRUxfo)) 
= £((V Vl ) T ) - B out (0i) - L(<h,vi). 
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So, the stress in the equation of vi can be rewritten as 



(mmx (Kf 2 + ]T f l KU j f k )) l 



i+j+k=2 

= ~((mm - h){Q f j2 + U ■ (0 l 1ZUcf> l ))) l - (mm x (f KU 1 fa)) l + L(&) 

= i((mm- il)(£((V Vl ) r ) -B ou M l ))) l - (mm x (JoKUxfa)) x + L(<f> u vi) 
= cri + (T 2 + L(^i, vi). 

Set a"3 = Di : (mmmm) f . Then the equation for vi can be rewritten as 

— — - — Avi + v • Vvi + Vpi = — — V • (<7i + a 2 + 03) 

- ^</oV?7i0 1 > 1 +V J &(0i,vi)+L(0i,vi). (7.13) 

In order to solve (17.12p -( j7.13p . we introduce the energy functional 

i?e 

£(t) = (4>i,4>i) + (01, C/i^i) + - (vi, vi). 

1-7 

Due to the choice of g, it is easy to see that 

(<h.,<h) + (<t>i, UkPi) > c((<k,<h) + <V0i, V^i)) . 
Notice that (V • (73, vi) = — (Di : (mmmm)j ,Di) < and by Lemma 18.81 and Lemma 18.61 

( - /C((V Vl ) T ) + BU<h), Uifr) - (V • (ai + a 2 ), Vl > 
< ( - /C((Vvi) T ) + -^(^(01) + Bcmtfai))) 

+ -<£ ou t(</>i),m ■ Vvi • m> + (U^^ui x (Vvi) T • m/ )> 

= ( - /C((Vvi) T ) + B in (<h.), -Aj^Bini^i) + BouMx))) + ^(BouM, -JU^/C (Di) + £(Di))> 

+ ( - Aj^BiM + fi out ((/)i)),/C((Vvif) +£((Wi) T )> 
= -{B in (fa),A] 1 B in (<f> 1 )) <0. 

Then by a simple energy estimate, we can deduce that 

±E(t)<C{l + E{t)). 

The estimate of the higher order derivative for (</>i,vi) can be obtained by introducing a 
similar energy functional. Once (fx, vi) is determined, we can get (/2, V2) and /3 by solving 
(|7.7p - (|7.8p in a similar way. Here we omit the details. □ 
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7.2. The remainder equations. We denote 

v = v + evi + £ 2 v 2 , f = fi + ef 2 +e 2 h, D = D + eD\ + e 2 D 2 , 
X T = (h + eh + e 2 h)TU £ h + hTU e (h + ef 2 + e 2 f 3 ) + f 2 TU e f 2 + f 3 T^^f 
, tn .f U e -U -eU 1 -e 2 U 2 t Us-lh-eUx , U £ - Up 

V £3 e 2 " 7 1 ^ g ^ 2 

f U e -U -eU 1 -e 2 U 2 -e 3 U 3 f U £ -U - eU r - e 2 U 2 
+ 70/1 n Jo H 13 /: 



e 



U e -U -eU l f U £ -U f ^ u £ -U -eU l f U £ -U f , 
H p / 2 H ~ 73 J + jYT ( ^ 70 H fl), ioxT = 1Z or V, 

Ll = " (it + V ° ' V/3 + Vl ' V(/2 + £h) + V2 ' V(/l + Sh + £2/s) ) 

+ ft • (x^ - m x (Vv ) T • m/3 - m x (V Vl ) T • m(/ 2 + ef 3 ) - m x (Vv 2 ) T • m(/i + ef 2 + e 2 / 3 )) , 

L 2 = "^r v • { D o : (mmmm/3) + Di : (mmmm(/ 2 + e/ 3 )) + D 2 : (mmmm(/i + e/ 2 + e 2 f 3 )) } 
1 — 7 

^-{V • (mm x X- R ) l + (Xy)i} - v 2 • Vvi - Vi • Vv 2 - ev 2 • Vv 2 . 

Then we can deduce the equations of (/jj, v^)(drop e for the simplicity): 

^jT + V • V/h + Vfi • V(/ + £/) + ^/oH/ /il 

= -U ■ (m x (Vv) T • mf R + m x (Vv R ) T • m(/ + ef) + e 3 m x (Vv R ) T • m/ H ) 

- ft • (f R KU £ f + f R K {Ue ~^ )k + JnU £ f R + e 2 f R 1ZU £ f R ) + L ls 

dv R _ „ o 7 

-— — + v R • Vv + v • Vvij + e°v R • Vv B + Vp - — Av fi 
at Re 

^-V • |d : (mm mm /'/,') 1 + D R : (mmmm(/o + e/))i + £ 3 D R ■ (mmram/ij)ij 



2Re 
1 - 

I-7//I 



i^V • (mm x (±f Q KH%f R + f R KU £ f + f R K ^ + f1ZU £ f R + e 2 f R RU £ f R )) i 

\foVH%f R + f R VUj+ f R V {U£ ~ £ U)f ° + fVUJa + e 2 f R VU £ f R ) ^ + L 2 . 



Re We 

Here D R = i(Vv fl + (Vv fi ) T ). We denote F R = F 1 H h F 6 with 

Fl = _ VK . V (/o + e7), 

F 2 = -n ■ (m x (Vv) T • mf R + f R 1W £ J+ /rJZ ^ ~ U ° )f ° + /TO^/h 



F 3 = -eTZ • (m x (Vv fl ) T • m/)) , F 4 = -e 3 TZ ■ (m x (Vv fi ) T • mf R ) , 
F 5 = -e 2 K ■ {fnWUfn) , F 6 = e 3 v R ■ Vf R , 
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and G R 



= Gi + 
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+ G% with 



~ ~ 1 — 7 ~ 

G\ = -vr • Vv - v • Vvr, G 2 = ~^T V - ( eD # : (mmmm/)i) , 

G 3 = • (mm x (D : (mmmm/^! + f R KUf + f R n {Ue ~^ )k + fRUfn))^ 

G 4 = l -^ 1 (-hVU £ f R + f R VU e J+ f R V {Ue ~ U)k + JvU £ f R ) , 
Re \e e 1 1 

1—7 

G 5 = ITB-^ ' ( D « : (mmmm/^i), G 6 = -e 3 v R • Vv R , 
2/te 

G 7 = -e 2 ^V • (mm x (f R KU £ f R )) v G 8 = -^^2(f R VUJ R ) r 
Then we rewrite the equations for (f R , v R ) as 

= -^-(mx(Vv R ) T -m/ )+F fi + L 1 , (7.14) 

IT " Re AVR ~ W V " (Dfi 1 ( mmmm /o)i) + Vp 

= • ( mm x {\kT^i%f R )) i + G fl + L 2 . (7.15) 



7.3. Some key estimates. In this subsection, we mainly present a control for a singular 
term in the error estimates. The proof is based on the lower bound inequality. Since we only 
have good lower bound for the part outside the Maier-Saupe space, we have to analyze the 
nonlinear interactions between the part inside the Maier-Saupe space and the part outside 
the Maier-Saupe space. Throughout this section, we will repeatedly use the notations from 
Section 5. Due to the assumption (|2.12p . we can construct a global coordinate transformation 
so that all results from Section 5 can be applied. 

Proposition 7.2. For any 5 > 0, there exists C = C(5) such that for any f £ i? 1 (Q x S 2 ) 
with Jg 2 /( x ) m)dm = 0, there holds 




We need the following lemmas. 



Lemma 7.2. It holds that 



Mil - M 22 = (a,2 



(A) - 2A 2 + Ai) 



2A 

2M 33 - Mu - M 22 = 3Co ' 



A A 4 - A\ 
At ' 
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Proof. Recall that M(x) = J g2 mm/dm. We get by direct calculations that 
=7: / sin 6>a dm + - / sm z 9a 2 dm = -Co — + Ca,2" 



2 ^ _ , 4 ^ " ^ u 2AI ' sa '^ 4A 

1 /" . 9 a , * 1 /" . 9 js , * , AA - M , (Aq - 2A 2 + At) 
M 22 =- / sin 2 6a drh - - / sin 2 6a 2 dm = -Co \\ 9 2 - 0,2- 



- 2 ^ _ 4 ^ _ su 2 ^ sa , z ^ 

-M33 = / cos 6laodm = Co , 

M12 =7 / sin 6>&2dm = 



4 J ' 4^o 

If., s. . . . A 2 - A 4 



1 /■ _ . , 

M13 =— / sin^ cos ^aidrh = Ca.i 

1 /" . - 

-^23 =2 / sin#cos#&idm = Cfe,i 



A 

Then the lemma follows. □ 
Lemma 7.3. There exists c > such that 

(n%f, f) > c{f\ f X ) + c(M kl , M kl - g £ * M u ). 

Proof. Noting that 

(H%f,f)=(H fo f,f)+a( f [ / (m.m') 2 /(x,m)/(x,m')dm'dmdx 

m/ (m • m') 2 g £ (x - x')/r(x, m)f R (x', m^dm'dx'dmdx) 
_ 2 J §2 / 

={U h f, f)+a I M(x) : (M(x) - g £ * M(x))dx, 
Jn 

then the lemma follows from Proposition 14.31 □ 
Lemma 7.4. We have 

\((fl-9e*fl),hf2) <c(i((/l- 5e * + \({h-9e*h\ /2> + </2,/2>). 

where the constant C depends on ||/3||l°° and ||V/3||x,oo. 



Proof. We write fi-g e *fi = (l - \{y/eD)) fi, that is x{0 = 1 - \A Hence, 

\{{H - 9e * fl), / 3 /2> =i<(l " x(Vi£»)) 2 /l, / 3 /2> 

=^((1 - x(^))/i, /s(l - x(V^D)f 2 ) 

+ i<(l - \WiD))h, [f 3 , X (V^D)]f 2 ). 
Then the lemma follows from the commutator estimate 

\\[f3,x(V^D)]f 2 \\ L 2 < C^||V/ 3 ||l-||/2||l2. 

By a scaling argument, it suffices to prove the commutator estimate with e = 1. Let Kj(pc) 
be the kernel associated with the Fourier multiplier i(djx)(D) (It is easy to show that Kj is 
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a Calderon-Zygmund kernel). We have 

3 



[f 3 , X (V^D)]f2 = V / #,-(x-y) / 1 5 J / 3 (rx+(l-T)y)dr/ 2 (y)dy, 

whose L 2 norm is bounded by || V/3||l°o ||j"2||l 2 ! see [2] for example. □ 
Now we are ready to prove Proposition [721 With the notations in section 5, we decompose 
/ as 

/ = ao(x, 6) + (afc(x, 6) cos kip + 6fc(x, 0) sin fc<£) . 
fe>i 

Then we can get 

-{fJd t {—))=-\ / <%(— - )(a (x,<9) + V] (a fe (x,0)cos/c^ + 6fc(x,^)sin/c( 1 5) ) dmdx 
£ Jo £ Jn Js 2 J° t>i 

^ C I I ~2 U + 7? + ft 2 + 7l + Pi + E( a ' + ^))dmdx 



InJsze k>3 



+ C / / (Co 2 + C 2 ,i + C 2 ,2 + C fe 2 1 + C b 2 2 ) dmdx 

+ - f f d t (-^r)((o(cos 2 8 - -r^) + sin0cos^(Ca,i cos(^ + Cb,i sin<^>) 
£ Jn Js 2 Jo v 

\ 2 

+ sin 2 #(Ca,2 cos 2<£ + C(, j2 sin 20) ) dmdx. 



As dt(l/fo) = —dtfo/fo an d ^t/o G Ker we may assume that 

11 > 

dt(~r) = 7(wi(x) cos </i> + t^2(x) sin 0) sin 6* cos 0. 
Jo Jo 



We have 



~[ [ d t (^)((o(cos 2 d - ^) + sin 6 cos 8 ((a t i cos ip + ( b i sin <p) 
£ Jn Js 2 Jo v M 

\ 2 

+ sin 2 8{Ca,2 cos 2<^> + (^2 sin 20) J dmdx 

2 f f 1 ~j4 
= _ / / )sin^cos(9(Ca,icos(^ + Cfe,isin^)(Co(cos 2 0- 

£ Jn Js 2 Jo A 

+ sin 2 #(Ca,2 cos 20 + Cb,2 sin 2<^))dmdx 

1 f f A 

+ - / / )(Co(cos 2 0- -^) + sin 2 ^(Ca, 2 cos 2(^ + 4,2 sin 2^)) 2 dmdx 

£ Jq Js 2 Jo A) 
2 /■ f „ , 1 , . A 5/ w _ o s A 2n 



/ )sin6»cos(9(Ca,icos^ + Cfe,isin^)(Co(cos 2 0- — 

£ Jn Js 2 Jo ^o 

+ sin 2 (?(Ca, 2 cos 2<£ + Cb,2 sin 2^))dmdx + - ( (Co + C,2 + C 2 , 2 )dx, (7.16) 

£ Jn 

where we have used the fact 

f 1 .... 2 

/ — (ioi(x) cos (/? + W2(x) sin <£) sin cos 9[ sin cos 8(C a l cos (/? + 0> l sine/))) dm = 0. 
Js 2 Jo 
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We get by Lemma 17.21 that 

2Co(x) - aN - a(2M 33 - M 22 - Mil - N ) 

= 2 Co (x) - 3aCo(x) ^ = MM _ M) Co(x). 

Note that the coefficient is positive in the front of Co( x ) if /o is a stable critical point. This 
implies that 

- f u; 1 (x)Ca,i(x)Co(x)dx< f 4(2Co(x)-aiV (x)) 2 + C(<5,||u; 1 || L oc)C Q 2 1 (x)dx 
e Jn Jn? 2 



+ j(2M 33 (x) - M n (x) - M 22 (x) - iV (x), C a ,i(x)wi(x)) 



On the other hand, we have 

Mij - Nij = A ki AijM kl - AkiAijige * M kl ) = A ki Aij(M k i - g £ * M H ), 

2 - 2 
Ca,l = —Mis = —A kll A vz M k , v , 
a a 

from which, Lemma 17.41 and Lemma 17.31 it follows that 
1 12 

-(Mij(x) - A^(x),Ca,i(x)wi(x)) = -(M kl - g £ * M kU -A^A^Av^v^w^)) 

< C(A, Wl )(-(M kl - g £ * M kl ,M kl ) + (M kh M kl )) 

<c(A, Wl )(- £ (n%f,f) + (fj)). 

Thus, it follows from Proposition 15.11 that 

- / Wl (x)Ca,i(x)Co(x)dx< / A(2Co(x)-aiVo(x)) 2 + C( ( 5,|h 1 || LO o)C 2 1 (x)dx 
eisi Jn? 2 

+ C (A,w 1 ){±(H%f,f} + (f,f}) 

< c(^(f nn%f,nn%f) + \{u%fj) + (/,/)). 

This gives the desired estimate for the term in ()7. 16p : 

2 /" /" 1 * -A 

- / dt(— ) sin^cos^Ca l cos c£>Co (cos 2 6 — — )dmdx. 

e Jn Jn 2 Jo A 

The other terms in (|7.16p can be treated similarly. We omit the details. □ 
The following lemma is used to control the other singular terms like -^v • V/r, Wfofii) in 
the error estimates. 



Lemma 7.5. We have 



||V/||| 2 < C{{n%Vf, Vf) + \{U e h f, /)) . (7.17) 



Proof. Let us first claim that 

\\f\\h<C{(njj) + (M[f],M[f])), (7.18) 
where M[/] = L 2 mm/dm. Due to the choice of g, we have 



K/(0i 2 < c(a - 5(ee))ime)i 2 + - e (i - 
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This implies that 

||V/||£ 2 < C«/ - g £ * V/, V/) + ±{f-g e *f, /)), 
which along with (|7.18p gives 

\\vf\\h <c{(n%vf, dj) + (vm[/], vm[/])) 
<c((n%vf,vf) + - e (H%j)), 

where we have used 

\^/ Q /> /) >« J J^ 2 f s2 ( m ' m ') 2 /( x > m ) (/( x > m - J ffe(x - x ')/( x '> m')dx') dm'dmdx 
=a<(l-«fe)*M[/],M[/]>. 

To complete the proof, it remains to prove the claim. We write / = /^ + / T with f T € KerC?f 
and / £ (KerC/^) -1 -. By Proposition 14.31 we have 

ll^lll^c^/,/)^^/,/). 

While, from the proof of Lemma 17.21 we know that 

||/ T |li 2 < C((M 13 ,M 13 ) + <M 23 ,M 23 » < C(M[/ T ],M[/ T ]> = C<M[/ T ], M[/ T ]>. 

This implies (I7TT8]) . □ 
In the nonlinear estimates, we will frequently use the following basic lemmas. 

Lemma 7.6. It holds that 

Mi (x)M 2 (x, m)M 3 (x, m) dmdx < C|| Mi 1 1 L 2 \ \M 2 1| #o, 2 1| M 3 1 \ L 2 , 

/ / Mi(x)M 2 (x,m)M 3 (x,m)dmdx< C||Mi|| h i||M 2 ||#o,i||M 3 || l2 , 
Jn J§ 2 

[ [ Mi(x)M 2 (x,m)M 3 (x,m)dmdx<C||VMi||^i||M 2 || i2 ||M 3 || i2 . 
Jn Js 2 

Proof. By Holder inequality and Sobolev embedding, we get 

/ / Mi(x)M 2 (x, m)M 3 (x,m)dmdx 
Jn Js 2 

< ||M 1 || L2 ||M 2 || L? o L2 ||M 3 (x,m)|| L2 < ||M 1 || i2 ||M 2 || i2i o ||M 3 || L2 



n Js 2 



<C||M 1 || L2 ||M 2 || H o, 2 ||M 3 || 



L 2 ■ 



The other two inequalities can be proved similarly. □ 
The following Bernstein type lemma is a direct consequence of Young's inequality. 

Lemma 7.7. Let k > be an integer and p > 2. Then it holds that 

||V fc Z4/|U P(nxS2) < Ce^^M^WfW^n^y 
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7.4. Error estimates. Let us first explain how to choose a suitable energy functional. It's 
helpful to look at the following toy model for (f R ,v R ): 

f + - £ A h n%f R = -K-{mx (Vvtff • m/o) , 



dt 



1 Avr + Vp = -^V • (mm x {±-f G KU%f R )) r 



dt Re ' " Re 
Compared with the homogeneous case, the new difficulty is caused by the singular term 
-V • (mm x (/o^'H/ /r)) 1 - To deal with it, it is natural to introduce the energy functional 

1 Re 

-(f R ,n%f R ) + y^(v r ,v r ), 

since we have the following important observation: 

(m x (Vv R ) T ■ mf ,KH%f R ) + ((mm x foR,H%f R )i,Vv R ) = 0. 

However, (/rj'H^/r) does not give a control for the part of f R inside the kernel. To have a 
control for the part inside the kernel, we need to introduce another functional (f R ,Aj^f R ) 
similar to the homogeneous case. So, the suitable energy functional for the toy model should 
be 

1 Re 
(fR,A£fR) + -(fR,n%f R ) + — ( YRt v R ). 

However, if we take ^(/fl,^/ /fl), a new singular term ^(f R ,d t (j^)f R ) will appear. Since 
there is no any decay in e for the part of f R inside the kernel, this term seems to have the 
order of 0(l/e)(Very singular!). Surprisingly, by analyzing the nonlinear interactions for 
(fn,dt(j^)f R ) and using the lower bound inequality, we find that it is bounded. 

In order to control the nonlinear terms, we also need to introduce a higher order analogous 
of the energy functional, whose choice is also very subtle. In all, our energy functional takes 
the form 

def/ 



<E e (t) ^(f R ,A£f R ) + -(f R ,U%f R ) + J^-(v*,v*> 



+C 1 e(Vf R ,A A lVf R ) + C 2 



+£ 3 (Af R ,Aj o 1 Af R ) + Q 



7 

e(Vf R ,n%Vf R ) + 



Re 
1^7 



e 2 {Vv R ,Vv R ) 



e 3 (Af R ,n%Af R ) + 



Re 
1^7 



e A ( Av R , Av R ) 



defl 



1 



= ^/*.«5b/*> + ^(n%fR,AH%f R ) + ^(Vv^Vv,) 

+ C l {Vf R ,H%Vf R ) + C 2 [(H%Vf R ,AH%Vf R ) +£ 2 ^— (V 2 v K , V 2 v fi ) 



+ e 2 (Af R ,H%Af R ) + C 3 e 2 (n%Af R ,An%Af R ) + e 



7 

4 7 
1-7 



VAvr,VAv r 



Here the constants C\, C 2 and C3 bigger than one will be determined later. 

Proposition 7.3. There exist c\ > and Eq > such that for any e £ (0,£o) an d t £ [0, T], 
there holds 

j£ e {t) + aUt) < C{1 + £ £ (t) + e^ZS?' 2 + e€ E {t?) + Ce 3 / 2 ^) 1 / 2 ^)- 
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where the constant C depends on ||/i||L°°(o,T;H 3 (r2x§ 2 ))(^ = 0,1,2,3), llv-tH^oc^T;// 3 ^))^ = 
0,1,2). 

Proof. From the definition of (£ E (t),$ E (t) and Lemma 17.31 it is easy to see that 

H/R|li a +C , i||e 1/2 / R |lHo.i + lk 3/2 /*fc 

+ ||vji||ia + C 2 \\£Vr\\ 2 h i + C 3 \\e 2 v R \\ 2 H 2 < C£ £ {t), 
llVvflUla + C 2 \\eVvr\\ 2 h1 + C 3 \\e 2 Vv R \\ 2 H2 < C&(t). 
And it is easy to show that 

||Ll||#0,2 + H-^HifS < C. 

These facts will be repeatedly used in the following calculations. For the simplicity of nota- 
tions, we denote A = Af and H £ = H 6 ^ in what follows. 

Step 1. 1? energy estimate 

Making L 2 (S1 x § 2 ) inner product to (|7.14p with Aj^fR, we get 

(^f^A^fR) + - e {n £ f R j R ) = -<v • Vf R ,A- l f R ) 

+ (m x (Vv fi ) r • m/ ,7^A-Vfl> + (Fr + L^A'^r). 
By Lemma 17.61 and Lemma 17.71 we have 

{F A ,A~ l f R ) < Ce l l 2 \\e 2 Vv R \\ H 4e l l 2 f R \\ H0 A\fR\\L* < Ce 1 / 2 ^ 2 ^), 
(F 5 ,A^f R ) < CV /4 ||/ R ||i 2 < Ce^£ £ (tf/ 2 , 
(F 6 ,A^f R ) < Ce\\f R \\ 2 L2 \\e 2 v R \\ H2 < Ce£ £ (t) 3 / 2 , 
and the other terms can be estimated as follows 

(v • Vf R ,A^f R ) + (F 2 ,A^f R ) < C\\f R \\ 2 L2 < C<£ e (t), 
(F 1 ,A~ 1 f R ) < C\\v R \\ L 2\\f R \\ 2 L2 < C(E £ (t), 
(F^A^fn) < Ce\\VvR\\ L 2\\fR\\ L 2 < C£ £ (t), 

So, we get 

(^A^fn) + ~ £ (n £ f R , f R ) < C(l + £ £ {t) + eV 2 £ £ (tf/ 2 ) + $&(*). (7.19) 
Make L 2 (tt x § 2 ) inner product to (j7TT4j) with U £ f R to obatin 

l(§- t fn,n £ f R ) + ^(nn £ f R ,f nn £ f R ) = ~(v • vf R ,n £ f R ) 

+i(m x (Vyr) t ■ mf ,nn £ f R ) + ~(Fr + L^Hefn). 
By Lemma 17.61 and Lemma 17.71 we have 

- e {F i + F^U £ f R ) < CJ^ftfR^Ho^VR^^nUefR^v < Ce 1 / 2 ^)^) 1 / 2 , 

- £ (F 5 ,H £ f R ) < Ce^WfRWlzWKHefRUi < Ce^ £ (t)Utf' 2 ■ 
Noting that 

/ V(/ + e/)dm = 0, / ATidm = 0, 
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we get by Poincare inequality that 

\{Fi,U E f R ) < C\\v R \\ L 4-KU £ f R \\ L 2 < C<Ee(t) + 6&(t), 

\{Ni,H s f R ) < C\\-KU £ f R \\ L 2 <C + 5de(t). 
We infer from Lemma 17.51 that 

-(v-vf R ,n £ f R )<c\\vf R \\ L2 \\-nn £ f R \\ L2 < -^=3 £ (t) + c<£ £ (t) + 5z £ (t). 

Here and what follows Go denotes a constant independent of 5. The other terms are estimated 
as follows 

\(F 2 ,U £ f R ) < C\\f R \\v\\±HH e f R \\v < C£ £ (t) + 6$ £ (t), 
\(Fz,Uef R ) < CellVv^ll^ll^J^II^ < Ced £ (t), 



Hence, we obtain 

-(^,H e f R ) + \(mn E f R ,KH e f R ) 

e ot ' e z N ' 

< C(l + <£ e (t)) + Ce l ' 2 Ut) 1,2 £e(t) + {S + e + -%)$&) 

+ i(m x (Vv fi ) T • m/o, 1lU e f R ). (7.20) 
Make L 2 (Q) inner product to ()7.15|) with v R to get 

2dt^ Rl Vi? ^ + Re^ VR ' VVjR ' > + ~2~Re~^ DjR : ( mmmm /o)i) , Vv fi ) 
= ^(mm x {- £ foKU e f R ), V Vji ) + (G fl + L 2 , v fl ). 

Obviously, (Gg, vr) = 0. We have by Lemma EH and Lemma 17171 that 

<G5, Vi? > < G^lle^/Rlko.ille^Rll^UVv^ll^ < Ce l l 2 £ £ (t)$ £ (t) l /\ 

<G 8 , v*} < G e 3 / 4 ||/^||| 2 ||v B |U 2 < Ge 3 / 4 e £ (t) 3 / 2 . 

While, (T) R : (mmmm/o)i, Vv_r) = \D R : (mmmm/o)i, > 0. The other terms are 
estimated as follows 

(G 1)Vii }<G||v^||| 2 <Gc£ £ (i), 

(G 2 ,v R ) < Ce\\Vv R f L2 < C<S e (t) + 5$ E (t), 

(G 3 ,v R ) < C\\f R \\ L 2\\Vv R \\ L2 < C£ £ (t) + 5d £ (t), 

<G 4 ,v R > < C^\\nn £ f R \\ L 2\\v R \\ L2 + C\\f R \\ L 2\\v R \\ L 2 < C£ £ (t)+5S £ (t). 
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Thus, we obtain 

^ VR ' Vfi > + ^ (VVR ' WR) 

< C(l + <£e(t) + e l/i £ e (t? /2 ) + Ce l l 2 £ e (t)$ e {t) l/2 + 

+ (mm x - £ hKU £ f R , Vv R ). (7.21) 

Step 2. 1 energy estimate 

Taking the derivative to (|7.14p with respect to Xj, then making L 2 (Q x § 2 ) inner product 
with eA~ 1 difn, we get 

e{^- t d t f R ,A~ l d t f R ) + (d t f R ,n £ d t f R ) 

= -{aifoRMeftonA^difn) - {di{^)f R ,dif R ) - e{di(v • V/hM" 1 ^) 

jo 

+ e<5 i (m x (Vv R ) r • m/o),^- 1 /^) + e(diF R + c^Li, .A"^/*). 
By Lemma 17.61 and Lemma 17.71 we get 

e(d i F 4 ,A~ 1 dif R ) < Ce\\£Vv R \\ H i\\e 1 ^f R \\ HO ,i\\^ 2 a i f R \\ H o,i < Ce$ e {t) l ' 2 £ e {t), 
e{d i F^A- 1 d i f R ) < Ce^ll/flll^He 1 ^!!^ < C e 5 / 4 £ £ (tf/ 2 , 
e{d i F^A- l d i f R ) < Ce^R^e^fdH^A^difR^ < Ced £ (t) 1/2 (B £ (t). 
It follows from Lemma 17.51 that 

e(di(m x (Vv R f • m/o),^- 1 ^) + s(d i F z ,A~ 1 d i f R ) 

^CIIeVvflllHill^/iillia < -%% £ (t) + 5$ £ {t)+C£ £ (t). 

VG2 

The other terms are estimated as follows 

(difoKHefaKA-idifR) + (di(^)f R ,dif R ) < CWfuhzWdifRh* < C€ £ {t) + 8$ £ {t), 

Jo 

- e<$(v ■Vf^^-'difn) +e(d t F 2 ,A~ 1 d i f R ) < Ce\\f R \\ 2 H o A < C£ £ (t), 
e(d t F 1 ,A^ 1 d i f R ) < Ce 1 /2||v fl || H1 || e V2a i / R || i . a < 8$ £ (t) + C£ £ (t). 
So, we get 

{j t dif R ,A- l dif R ) + - £ {u £ dif R ,d t f R ) 

< C(l + <E e (i) + e^^Sf 2 ) + Ce$ £ (t) 1/2 £e(t) + (6 + (7.22) 

V^2 

Taking the derivative to (|7.14p with respect to Xj, then making L 2 (f2 x § 2 ) inner product 
with eUdifn, we get 

£(^- t dif R ,H £ dif R ) + (f nu £ d i f R ,ii'H £ d i f R ) 

= -(difoRLUfR + fRRS i (Uofo),'R,'H e a i f R ) - e<^(v • Vf R ),H E dif R ) 
+ e($(m x (Vv fl ) T • mf ),nn e dif R ) + e<SiF fl + d^U^f R ) . 
The first term on the right hand side is bounded by 

H/flllialiroWfllli,, < ce £ (t) + 
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By Lemma 17.61 and Lemma 17.71 we get 

e{d i F A ,H e d i f R ) < Ce^e^RMe^fRUo^TZHedifR^ < Ce^ 2 <a £ (t)d £ (t) l/2 , 
e( diF^KUedifR) < Ce^H^/ifJkoall/ifJI^H^^/Hll^ < Ce 7 l 4 £ £ {t)$ £ {t) l '\ 
e{d i F^n e d l f R ) < Ce l l 2 \\^-v R \\ H , || £ 3 / 2 /^||ho, 2 ||^ £ ^|| l2 < Ce^ 2 <a e (t)de(t) l/2 . 

and by Lemma 17.51 

e{di{v-Vf R ),U E dif R ) 

= e(div ■ Vf R , ^p-) - e(v • V(^)dif R , dJ R ) - e(v • Vf R , 8^8^)) 
Jo Jo 

< C^WVfnWvWaifRWv < C£ £ (t) + 5$ £ (t). 
The other terms are estimated as follows 

e(a i F 1 ,n e a i f R ) < aiieVvflii^ii^^/^n^ < cg^t) + 

e{d i F 2 ,H e d i f R ) < C7e 1 / 3 ||e3/ Ji || H o,i||^ e ^/fl||i a < C<£ e (t) + 6$ £ (t), 

s(diF 3 ,H £ dif R ) < CWe^RW^WKnedifnW^ < C£ £ (t) + 6$ £ (t), 

e(m x (V Vi? ) T • mdifanHedifn) < C|| e v R || H i \\nn £ dif R \\ L 2 < C£ £ (t) + 53 £ (t). 

Thus, we obtain 

e(j t dif R ,u £ d t f R ) + (f nn £ d t f R ,nn £ d t f R ) 
< c(i + e E (i)) + c^ 2 ^{m e {t) 112 + 

+ i(m x (VdiV R ) T • xxif ,nH £ dif R ). (7.23) 
Making L 2 (£7) inner product to ([7. 15ft with s 2 dfvji, we get 

~2dt^ 9{VR,diVR ^ + ^ e2 ( V<9 * Vi? ' Va ' V/? ) + ~2Re £2 ^ di ^ R : ( mmmm /o)i), Vd{v R ) 
= ^e{mm x d i (f Q nH £ f R ),Vd i v R )+e 2 {d i G R + d i L 2 ,d i w R ). 

First of all, we know that ( s (diD R : (mmmm/o)i), V^vr) > and 

e 2 ((U R : ^(mmmm/o)i),Va iVii > < C\\ev R \\ H i ||eV^v fl || £a < C<£ £ (*) + 5&(t). 

By Lemma 17.61 and Lemma 17.71 we get 

e^Gfc.ftv*) < Ce^We^fuWHO^RWH^VvjiW^ < Ce?<£ e (t)Ut) 1/2 , 
e^^Ce,^) < Ce\\e 2 -v R f H2 \\Vv R \\ L 2 < Ce£ £ {t)$ £ {t) 1 ' 2 , 
e 2 (d l G 7 ,d l ^ R ) < Ce 3 / 4 ||e 1 /2 /ii || H o,i||/R|| ia ||e 2 v ii || fla < Ce^^tf' 2 , 



e 2 



(diGs,diV R ) < Ce\\f R \\ L 2\\e 1/2 f R \\ H o,i\\e 2 v R \\ H 2 < 
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And the other terms are estimates as 

e 2 {diG u diV R ) < C\\ev R \\ 2 H i < C£ £ (t), 

e 2 (d l G 2 ,d i v R ) < C|| e 2 v R ||^ 2 ||eAv R || L 2 < Ce<£ e (t) + 8$ e (t), , 
e 2 {d i G^d i w R ) < Ce^ 2 \\e l l 2 f R \\ Ha ,4eVd^ R \\ L 2 < C£ £ (t) + 6$ e (t), 
e 2 (d l G,,d i v R ) < C\\f R \\ L2 \\eAv R \\ 2 L2 < C£ e (t) + 6$ e (t). 

So, we get 

£ 2 (^diV R ,diV R ) + e 2 r^^{VdiV R ,VdiV R ) 

< C(l + £ £ (t) + e 3/i £ £ (tf/ 2 ) + Ce^ £ (t)d £ (t) 1/2 + *&(<) 

+ (mm x ^foKU £ dif R , Vd l v R ). (7.24) 

Step 3. H 2 energy estimate 

Taking A to (|7.14j) . then making L 2 (Q x S 2 ) inner product with e 3 A~ 1 Af R , we get 

e^A/a^A/fl) + e 2 (Af R ,U £ Af R ) 

= e 2 {f R KAU fo + Af -RU £ f R + 2d i f KU £ d l f R + 2$/ fl 72W 3i/o, A//?> 
- e 3 (A(v • V/aM^A/fl) + e 3 (A(m x (Vv^f • m/o^T^A/fl) 
+ e 3 (AF R + AL 1 ,^- 1 A/ i? >. 

The first term on the right hand side is bounded by 

\\e?/ 2 KA- 1 Af rHIs + We^difRWl* +s\\f R \\ 2 L2 < C£ £ {t). 

By Lemma 17.61 and Lemma 17.71 we have 

e 3 <AF 4> ^- 1 A/ H > < Ce||e 2 Vv i? ||^||6 3 / 2 A/ i? || L2 || e 3 / 2 / i? || H0 , 2 < Ce^S^Sf'^ 

e\AF^A- l Af R ) < Ce">l A \\f R \\ L 4f R \\ 2 m2 < Ce^^tf' 2 , 

e s (AF &1 A- l Af R ) < a e || £ 2 Vv^|| J , 2 || £ 3 / 2 A/ K || i2 || e 3 / 2 / Ji || J ,Q, 2 < Ce<E e {t)^{t) 1 ' 2 . 

And by Lemma 1731 the term e 3 (A(Fi + F 2 + F%), A~ l Af R } is bounded by 

e^Wevr \\ H 2 ||e 3 / 2 Af R \\ L 2 + \\e z ' 2 f R \\ 2 m , + \\e 2 Vv R \\ H , \\eAf R \\ L 2 



<c<*s) + {^= + t)Ut)- 



So, we get 



{-Af Rl A- l Af R ) + -(H £ Af R ,Af R ) 

< C (l + <S E (i) +e 5 ' A £ £ {tfl 2 ) +Ced £ (t) 1/2 £s(t) + (7.25) 

V^3 
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Taking A to (|7.14p . then making L 2 (f2 x § 2 ) inner product with e 3/ HAf R , we get 

£ 3 {^- t dif R ,n E Af R ) + e 2 (f nn £ Af R ,nn £ Af R ) 

= e 2 (f R KAU f + Af KUJ R + 2d i f Q VU e d i f R + 2dJ R TZU dJo, KH E Af R ) 
+ e 3 (A(m x (Vv fl ) T • mf ),1lU £ Af R ) + e 3 (AF R + AL u U £ dif R ). 

The first term on the right hand side is bounded by 

|| £ ^H £ A/ Ji || L 2( e 1 / 2 || e 1 /2a. /K || i2+£ || /R || L2 ) < 6$ E (t) + C£ e (t)). 

By Lemma I7.6I and Lemma I7.7I we get 

e 3 (AF 4 ,n £ Af R ) < Ce 3 / 2 ||e 3 Vv fl || fl3 ||e 3 / 2 /ji||H0, a ||eWW e A/ ii || La < Ce 3 / 2 3 £ (t)<Se(t) 1/2 , 
e 3 (AF 5 ,KH £ Af R ) < Ce 7 / 4 \\f R \\ L 2 \\e 3 / 2 f R \\ H o,2 \\eK% e Af R \\ L 2 < Ce 7 l A £ £ (t)$ £ (t) 1 ' 2 , 

and for Fq, we have 

e 3 (AF 6 ,U £ Af R ) = £ 6 (Av, • Vf R , ^> + 2s 6 (d l v R • V^/r, ^> 

Jo Jo 

+ £ 6 {v fl • Vf R , A(U £ Af R )) -^-(v R - V(1)A/ R , Af R ) 

z Jo 

<Ce{\\e 2 Vv R \\ H 2\\e 3 / 2 Vf R \\ H o,i\\e 3 / 2 Af R \\ L 2 

+ e 7 / 4 ||v fl || L2 ||e 1 / 2 V/ i? || L2 || e 3 / 2 A/ i? || L2 + || e 2 v fl || H2 || £ 3 / 2 A/ fi || 2 2 ) 
<Ced £ (t) 1/2 <B £ (t) + Ce(B £ (t) 3 / 2 . 

And the term e 3 (A(Fi + F 2 + F 3 ),H e Af R ) is bounded by 

II^VRll^lle^A/^ll^+e^ii^^ii^ii^^^n^ 
+e||e 2 V Vi? |^ 2 || £ ^H £ A/^|| L2 < C£ £ (t) + (e + <*)&(*), 

and by Lemma 17.71 

e 3 (A(v • V f R ),1i £ Af R ) = e 3 (Av • V/*, ^> + 2e 3 <d l v • Vft/*, ^> 

7o 7o 

+ e 3 (v • V/ flj A(Z4A/r)> - y <v • V (^) A /«> A /«> 

< CII^v/rIU.II^A^II 2 , +C||, 3 / 2 A/ r ||| 2 < <?<£,(*)• 

So, we get 

e 3 (^f R , n £ Af R ) + e 2 (nn £ Af R , f nn £ Af R ) 

< C(l + C e (t) + e£ £ (t) 3/2 ) + C £ (E £ (i)S £ (i) 1/2 + C(<J + £ 3/2 <£ £ (*) 1/2 )S £ (t) 

+ e 3 (m x (VAv fi ) T • mf ,nn £ Af R )). (7.26) 
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Making L 2 (S7) inner product to (|7.15p with £ 4 A 2 vr, we get 
^--(Av R ,Av R ) + ■Xe 4 (VAvj i ,VAvj i > 



1-7.4 



2i?e 



e 4 A D R : (mmmm/o)i), VAvfl 



+ L_L £ i( mm x A(f TZn £ f R ) , VAv R ) + e 4 (AG R + AL 2 , A V/? >. 

Again, ((ADjj : (mmmm/o)i) , VAv^J > 0. The other part of the first term is bounded by 

H^Vflll^H^VAvflllia < C<£ £ (t) + 6dS). 

By Lemma 17.61 and Lemma 17.71 we get 

e 4 (AG 5 ,Av fi ) < C f e 3 / 2 ||e 3 / 2 / R ||H0.3||e 2 VAvj i || £2 ||e 2 VvB|| fla < C^ 2 ^) 1 / 2 ^*), 
£ A {AG 6 ,Av R ) < Cs\\e 2 v R \\ 2 H2 \\s 2 VAv R \\ L2 < Ce<S E (t)$e(t) 1/2 , 
e 4 <AG 7 ,Av R > < Ce 7 / 4 ||/R|| La ||e 3 / 2 / R ||flo, 2 ||e 2 VAvj i || La < Ce 7 ^£ £ (t)$ £ (t) 1 ' 2 , 
e 4 (AG s ,Av R ) < Ce^We^fuW^WfRh^VAvnWv < s<£ £ (t)Ut) 1/2 ■ 

And the term e 4 (A(Gi + G 2 + G 3 + G 4 ), Av R ) is bounded by 

\\e 2 v R \\ 2 H2 + e\\e 2 Vv R \\ 2 H2 + e 1/2 ||e 3/2 /fl||ff°.3 ||e 2 VAv i? || L2 < C<£ £ {t) + (5 + e)$ e (t). 
Then we get 

< C(l + <E £ (t)) + (5 + e)&(i) + CeZS^S) 1 ' 2 

+ Ce 3 / 2 ^) 1 /^*) + (mm x -J Q n% £ Af R) VAv R ). (7.27) 

Step 4. The closing of the energy estimates 
Noting that 

-(m x (Vv*) T • mf ,KH £ f R ) + i((mm x / KW £ /r)i,Vv b ) = 0, 

and then summing up (I7.19j) - (l7.27p . and taking Ci big enough, and then C%, C3 big enough, 
and finally taking 5 small enough, we infer that there exist e > and c\ > such that for 
any e E (0, £0), there holds 

(^/R.^ _1 /fl> + \{^- t fR,K £ f R ) + r^i^VR, V R ) 

+C 1 £(-V/ i? ,^ 1 V/ fl > + C 2 e{-Vf R ,H £ Vf R ) + C 2 e 2 r ^(-Vv R , V Vi? > 

+£ 3 (^A/ R ,^ 1 A/ i? > +C 3 £ 3 <^A/ i? ,^ e A/ i? > + C 3 £ 4 r ^(^Av K ,Av i? ,> + ci&(t) 
< C(l + <E e (t) + ^(^(i) 3 / 2 + £(£ e (t) 2 ) + Ce^tS) 1 /^). 
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Now Proposition 17.21 implies that 

~(f R ,n £ f R ) < (^,H £ f R ) + 5Ut) + C£ £ (t), 

~(Af R ,H e Af R ) < (?^,H £ Af R ) + 6Ut) + C£e(t), 
and we have the trivial inequality 

~(f,^f)<(p,A^f) + C\\nA^f\\h. 

Thus, we can deduce that 

j£ E {t) + aUt) < C{1 + <E e (t) + e^Zeitf' 2 + e£ £ {tf) + Ce 3 / 2 ^) 1 / 2 ^*). 

This completes the proof of Proposition 17.31 □ 
Now we are ready to prove Theorem 12.31 Given the initial data (/g, Vg), we can show by 
the energy method [23] that there exists T £ > and a unique solution (/ e (x, m, i), v e (x, t)) 
on [0,T e ] to dUHD-dUl]) such that 

f e (t) - 1 G C([0,T £ ];H\n x S 2 )), v £ (t) G C([0,T £ ]; H 2 (n)) D L 2 (0,T £ ; H 3 (n)). 

While, Proposition 17.31 tells us that 

j&S) + dff e (t) < C(l + + e 1 / 4 ^) 3 / 2 + 6£ e (t) 2 ) + cV 5 / 2 ^) 1 / 2 ^), 

for any t G [0,T £ ]. Due to the assumptions of Theorem 12.31 we know that (£ e (0) < C. Thus, 
there exist sq > depending on T such that for any e G (0, Eq) and t G [0, min(T, T e )], there 
holds 

CeW+Cl / S £ (s)ds< C. 

This in turn implies T e > T by a continuous argument. Then Theorem 12.31 follows. □ 



8. The dissipation of the Ericksen-Leslie energy 
Recall that the Ericksen-Leslie equation has the following energy law 

-U [ Re J v| 2 dx + ^) 
d*Wn2(l-7)' / 

= / I—?— |Vv| 2 + (Q 1 + ^ 2 -)|D:nn| 2 + a 4 D:D 

Jn v 1 - 7 7i 

-Y 2 1 \ 

+(a 5 + a 6 - — )|D • n| 2 + — In x hi 2 dx. (8.1) 

7i 7i ' 

Because the relations between six Leslie coefficients are unclear in Physics, whether the 
energy is dissipated remains open. In |12j . Lin and Liu present some constrains on the Leslie 
coefficients to ensure that the energy is dissipated. We will show that the energy (|8.ip is 
dissipated for the Ericksen-Leslie equation derived from the Doi-Onsager equation. More 
precisely, 
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Theorem 8.1. If the Leslie coefficients are determined by \2. 6\) and (2.1), then there holds 

2 2 

(ai + — )|D : mil 2 + a 4 D : D + (a 5 + a 6 - — ) ID • n| 2 > 

7i 7i 

for any symmetric matrix D and n 6 S 2 . 

Remark 8.1. Recall that 71 = S2/A. By taking u = u' in \8. 3\) . we see that A > 0, thus 
71 > 0. 

Throughout this section, we denote by /o = fa^n a critical point of A[f]. 

8.1. Some useful identities. Recall that S2 = (P2( m-n ))/ an( i ^ = (-f4( m • n ))/ > where 
Pk(x) is the fc-th Legendre polynomial. We define 

M*- 4 ^ = (mm mm) * . 

Lemma 8.1. It holds that 

M( 2 )=5 2 nn+i^I, 

S4 2S2 1 

+ naTtpSfa + npn^da^) + (— - — + — ) (8 a p8 7IJi + eVy^ + 8 ail 8^). 

The lemma is a direct consequence of Lemma 16.51 Especially, the lemma implies that 
Lemma 8.2. For any symmetric matrix D ; there hold 

M (2) • D =5 2 n(D • n), D M (2) = 5 2 (D • n)n; 

M( 4 ) : D =<S 4 nn(D : nn) + 2 ^ ~ S ^ (( D . n ) n + n ( D . n )) 

+ 2(— + — )D + — -ID : nn . 

V 35 21 15 ; 7 V ' 

Lemma 8.3. For any symmetric constant matrix D, there holds 

(K ■ (m x D • m/ ) , f) = \t> : / (mm - \l)K ■ (/o^/)d: 

* Js 2 -J 



m. 



Proof. It is easy to show that for any vector field v defined on § 2 , 

((mm \)TZ ■ (/v))i = ((m x v)m + m(m x v))/. 

3 

Applying it with v = IZg and v = m x (k • m), we deduce that 

r 1 

/ (mm 1)7?. • (fTZg)dm = (m x IZgm + mm x IZg) f. 

J§ 2 3 

Thus, we have 

(K ■ (m x D • m/ ), /> = D : (m(m x Kf)) fo 
= Id : ((m(m x Kf)) fo + <(m x Kf)m) fo ) 

= l -T> : f (mm - \l)K ■ (/ ft/)dm. 
The lemma follows. □ 
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Lemma 8.4. For any antisymmetric constant matrix fi, we have 

K ■ (m x (O • m)/ ) - (n x (O • n)) • Kf = 0, 

Proof. The lemma is a direct consequence of the following identities 

K ■ (m x (ft ■ m)) = K i (e ijk rn j tt k im l ) = (I - 3mm) : CI = 0, 
(m x (ft ■ m)) • IZfo = (m x (n • m)) • (m x n)/g 
= (n x (O • n)) • (m x n)/£ = (n X (O • n)) • Kf . 

The proof is finished. □ 

8.2. Projection operator and properties. We denote by Pi n the projection operator from 
Vq(E> 2 ) to Kert?y , and denote by P ou t the projection operator from Vo(S 2 ) to (Ker^)^. Since 

KerGf is orthogonal to (KeiGf ) ± under the inner product (-,Aj (•)}, we have 

(f,Aj*f) = (Pin/^^Pin/) + (Pout/^^Pout/). 

For any constant matrix k, we define 

K(k) = P in [K ■ (m x (as • m)/ )] , C(k) = P out [ft • (m x (as • m)/ )] . 
Lemma 8.5. It holds that 

£(k) = (n x (AD • n - Q ■ n)) ■ Kf . 
Here D = \{k + k t ), ft = ±(k t - as). 
Proof. By Theorem l4.lt we may assume that 

P b [K.(mx(K-m)/ )]=w-K/o 
for some vector w with w_Ln. Thus for all • 7£/o € Ker<5y , 

(ft ■ (m x as ■ m/ ),^7 o 1 (e • ft/ )> = (w • IZfaAj^Q ■ Ufa)). 
First we claim that 

(wKf ,A] \&-Kf )) = &-{nx (S 2 D-n-^ft-n)). (8.2) 

Let u and u' be any vectors. By Proposition 14.41 we may write 

A~J*(u' ■ ft/ ) = (u[ sine/)- u' 2 cos <j))g {9). 
Then we get by a direct computation that 

= / 2n sin 9 cos 6fo(ui sin (j> — U2 cos (/>){u'i sine/) — u' 2 cos <fi)g(8)dm 

= J(ux n)(u' x n) f f ^-g(6)dm = -l(u - (u ■ n)n) • u'. 

2 J S 2 d6 A 

Therefore, w = n x (AD • n — • n). 



(8.3) 
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Now, we prove f|8.2j) . By Lemma IH731 Lemma [83] and (|8.3p . we have 
(H-(mx K -mf ),A-\&-Kf )) 

= (ft • (m x D • m/ ),^- 1 (e • ft/ )> - (ft • (m x ft • mf ),A-\® • ft/ )> 
= ~((mm : D), • ft/ ) - <(n x (O • n)) • ft/o,^© • Ufa)) 

= (0 • (m x (D • m))) /o - ^0 ■ (n x (O • n)) 

A 

= S 2 e • (n x (D • n)) - ^0 • (n x (ft • n)). 

A 

The claim follows. □ 
Lemma 8.6. C(fl) = for any antisymmetric matrix ft. 

Proof This is equivalent to prove JC(ft) = ft ■ (m x (ft ■ m)/o), which is a consequence of 
Lemma 18.41 and Lemma 18.51 □ 



Lemma 8.7. For any symmetric matrix D ; there holds 

f £(DM 7 ;£(D)dm=(^±^i-A5 2 )|D.n| 2 

_ ,54 _ 252 U {xs _ s )( nn)2 _ 

v 35 21 15 ; 1 n ; 

Proof. Applying Lemma [8731 with / = Aj^(TZ ■ (m x D • m/o)) and Lemma [8721 we get 

(ft • (m x (D • m)f ),A]*K • (m x (D • m)/ )) 

If 1 
= D : / (mm )ft • (m x (D • m)/o)dm 

2 y§2 3 

= — — D : (2D : (mmmm) f — D • (mm) j — (mm) f ■ D) 
= -D : (S 4 nn(D : nn) + 2 ^ 2 ~ ^ ((D • n)n + n(D • n)) 

+ 2(|-f + > + ^.(D:nn,) +S2 (D.n)^ 

35 2 + 45 4 , . 2 ( Si 2S 2 1 x_ _ cm \2 

= (D • n) — 2 D : D — 64 (D : nn) , 

7 v ' v 35 21 15 ; v ; 

which along with Lemma 18.51 gives 

{C(Ii),Af £(!>)) = <P out (ft • (m x (D • m)/ )), ^Pout {n • (m x (D • m)/ ))> 
= (ft ■ (m x (D • m)f ),A^K • (m x (D • m)/ )) - AS 2 |n x (D : n)| 2 

= (™*±±S± - A5 2 )(D • n) 2 - 2(| - ^ + 1)D : D + (XS 2 - 5 4 )(D : nn) 2 . 
The proof is finished. □ 
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8.3. Proof of Theorem 18.11 and application. Let us first prove Theorem 18.11 By (|2.6j) - 
(12.71) and Lemma 18.71 we find that 



2 2 

(at + ^)|D : nn| 2 + q 4 D : D + ( 73 - ^)|D • n| 2 

7i 7i 

= (-§ + — )|D : n„|* + (-| - ^ + A )D : D + (5*+^ _ i)| D . n |2 
2 71 35 21 15 7 71 

^^^-|)(D.n) 2 -2(|-§ + l)D : D + (- S4 + i)(D:nn) 2 

+ ^^(D • n) 2 + |(D : nn) 2 + (| - ^ + jj)D : D (8.4) 



>(^^-^)(D-n) 2 -2(|-^ + l)D 2 + (-S 4 + i)(D:nn) 2 >0, 
7 71 35 21 15 71 

since all the coefficients in the line (|8.4p are positive. Indeed, we have 



c /-'■/o^ \2 i\\ 3A 2 - A) 
S 2 = (-(3(m-n) - l)) /o - 



l 2 v v ' 2A 
5 4 = <^(35(m • n) 4 - 30(m • n) 2 + 3) 



l /or/ , 4 2 , nxX 35A4 - 30A 2 + 3A c 



A 8A 

(A 8 - 2A 6 + A A ) > 0. 



8^0(277) 



Hence, 



S 2 -S 4 = J r (6A 2 - 5A 4 - A ) = — — ^ — (Aq - 2A 4 + A 2 ) > 0. 

This complete the proof of Theorem I8.ll □ 

As a byproduct, we get the following dissipation law, which has been used in the existence 
of the Hilbert expansion. 

Lemma 8.8. For any matrix k, there holds 

((mm-h)C( K )) 1 : k< 0. 

Proof. Lemma I8.3I implies that 

A fo (D : (mm - -I)) = -2K • (m x D • m/ ) = -2(/C(D) + £(D)). 
3 

Here D = ^(k + k t ). From Lemma \8.6\ we know that C(k) = C(D). Hence, 

1 f 1 
((mm- -I)C(k))-. k = / (mm 1) : D£(D)dm 

3 y§2 3 

= -2<^ ) 1 £(D), /C(D) +£(D)> = -2<^ o 1 £(D), £(D)> < 0. 
The proof is finished. □ 
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